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The Effect of Variation of Mass on the 
Dynamic Stability of Jet-Propelled Missiles’ 


M. V. BARTONT 
The Unwersity of Texas 


ABSTRACT 


The variation of mass of a jet-propelled missile, caused by ex- 
pending fuel, produces jet damping forces that, in addition to the 
variation of the mass itself, affect the dynamic stability and, 
hence, the flight behavior of the missile. The equations of mo- 
tion in the vertical plane for an airjet-type missile and a rocket- 
type missile having a constant rate of loss of mass are developed, 
and the characteristics of the jet damping forces are indicated. 
The equations of motion are solved for the simplified case of a 
symmetrical rocket missile in free flight, and an evaluation is 
made to determine the effect of a gust disturbance on the subse- 
quent behavior of the missile. These results are compared with 
the behavior of a constant mass missile. It 
disturbance of the variable mass missile damps out more rapidly 
than it does for the constant mass missile, indicating that the 
The oscillation frequency 


is shown that the 


variable mass missile is more stable. 
for the variable mass missile increases with time and is greater 
than the oscillation frequency for the constant mass missile 


INTRODUCTION 
AS A JET-PROPELLED MISSILE 


carries its own oxygen or an airjet that takes 
expends fuel, its mass 


such as a rocket that 


oxygen from the atmosphere 
decreases. The effect of the loss of mass has been con- 
sidered in analyses of the performance of high-altitude 
+fockets and in the problem of escape from the earth by 
many investigators, among which are Malina and 
Smith,' Tsien and Malina,? Malina and Summerfield,* 
and Seifert, Mills, and Summerfield.‘ In these cases 
the dynamic stability of the missile has not been of par- 
ticular concern. For rocket- or airjet-type missiles 
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that are to be flown either guided or unguided mainly 
within the atmosphere, the problem of stability and 
control is an important consideration. Many of these 
missiles have aerodynamic surfaces such as wings and 
tails whose function is to control and stabilize the mis- 
sile in flight. It has been determined that the variation 
in mass causes so-called jet damping forces that should 
be considered along with the aerodynamic forces in 
determining the stability of missiles. The effect of 
these forces caused by the variation in mass will be con- 
sidered herein. 

There are a number of ways of expressing the equa- 
tions of motion for a system of variable mass. One of 
these methods involves a study of the change of mo- 
mentum in an infinitesimal increment of time, and an- 
other is a study of the effective forces acting at any 
instant. Considering the first of these methods, it is 
possible to interpret Newton's equation as the equality 
between the change of momentum and the impulse 


for an interval dt, or 


d(mV) = Fdt 


The change in momentum should account for the change 
of momentum of the initial system, as well as the 
change of momentum of the particles leaving the system. 
An understanding of the kinematics of the motion is 
necessary for successful application, since relative 
velocities between detaching particles and the parent 


mass must be considered. 


The second method, which is the one used here, is 
particularly suited to systems for which the mass flow 
is steady. This method depends on the equality be- 
tween the “‘effective forces’’ of the particles of a sys- 
tem and the external forces. The mass of a particle 
times its acceleration is called the effective force. The 
physical significance of the jet damping forces is clearly 


indicated by this method. 
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Fic. 1. Missile. 








Fic. 2. Rocket jet. 





Entrance Plane- 











Fic. 3. Airjet flow. 


{QUATIONS OF MOTION 


A schematic diagram of a missile is shown in Fig. 1. 
Two propulsive arrangements are to be considered: 
the rocket unit shown in Fig. 2 and the airjet unit 
shown in Fig. 3. The following assumptions will be 
made. 


(1) Only motion in a vertical plane is considered. 


(2) In the case of the rocket, the system investi- 
gated at any instant consists of the missile body or 
frame, the jet inside the body from the exit plane to the 
fuel injection point, and the fuel. The airjet consists 
of the foregoing components plus the airflow inside the 
missile from the entrance plane at the nose to the fuel 
injection point. 


(3) The external forces are lift, drag, and weight. 
The aerodynamic forces are caused by pressures that 
are measured from an atmospheric datum. 


(4) The reservoir of fuel is located symmetrically 
about the longitudinal axis of the missile. The fuel 
has no velocity relative to the missile frame. 


(5) The mass flow is symmetrically located about 
the longitudinal axis of the missile. The mass flow is 
steady—that is, the mass of material passing any cross 
section at any time is a constant. 


(6) The missile is roll-stabilized, and cross-coupling 
terms between a symmetrical motion (roll, yaw, an 
sideslip) and a symmetrical motion (pitch, longitudiny 
displacement, and transverse displacement) are no 
considered. 


Consider the missile shown in Fig. 1, with the X ay 

Z axes fixed relative to the body and with the origi 
at C. If the missile is rotating with an angular velocit, 
Q, then the acceleration of any particle A of the rigi 
frame is the vector sum of the absolute acceleratig, 
Ve of the reference point C plus the tangential acceler. 
ation Or and the normal acceleration Q?r of A relatiy, 
to C. Since a particle B of the interior jet has a ye 
locity v and an acceleration @ relative to the rotating 
axes XZ, its absolute acceleration will have the accel. 
eration # in the negative X direction and the Coriolis 
acceleration 2Qv in the Z direction, in addition to the 
accelerations of the type indicated for a particle of the 
frame. 

The effective force for a particle is the product of its 
mass and the absolute acceleration. The effective 
forces for the particles of the frame and jet consist of 
the effective forces for both the frame and jet, as if 
these particles composed a rigid body, and the effective 
forces of the jet caused by the motion of the particles 
of the jet relative to the rigid body. Equating the 
sum of the effective forces and moment to the external 
forces and moment, we have 


Vox dom, oa OF mers - OS mr - oe = i. 
Feed — OLmey — OL mee + 

20¥m~e = F, ) 
OS + i na = Me 


where 
rx, tz = X and Z components of distance between 
any particle 7 and C 
b 


Z. 


summation over all particles of body 1 
cluding frame, jet, and fuel at any i 


stant 

J 

> > = summation over all particles of jet at any 
instant 

V = dV /dt, etc. 


Choosing point C as the instantaneous center 0! 
mass, then 
b b 
yomirx ss om irz = 0 
The instantaneous mass moment of inertia about ¢ 
is 
b 
Ym =I 


The components of V¢ in directions of the rotating 
axes XY and Z are 
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DYNAMIC STABILITY OF 
Vy + V2Q in the X direction | (2 
Vz — VxQ in the Z direction { " 


Using C as the center of mass and substituting the 
values of the components for Vc into Eq. (1) gives 


mVx + mV,Q- Smo = Fy 
j 
mVz — MVxQ + 20> my = F, 
j 
10 + 20>-mipv = Me 


where m = total mass of missile including the jet at 
any instant. In the first of these equations, the term 
involving the jet acceleration is the jet thrust, while 
the jet terms in the last two equations are the jet 
damping force in the Z direction and the jet damping 
pitching moment, respectively. 

In the study of stability and control of aircraft, it is 
usually convenient to consider only the small changes 
of attitude, velocities, and accelerations from a steady- 
state flight condition. These changes are called per- 
turbations. The quantities are expressed in terms of a 
constant part plus a small change. Thus, for ex- 
ample, the velocity in the X direction is expressed as 
the sum of a constant velocity and a change in veloc- 


ity. 


U,+u 


Vx = 
similarly, 
Q0=Aa+ 9 


The lower case letter usually denotes the change in the 
quantity. This follows the notation used in Durand’s 
Aerodynamic Theory.2 The external aerodynamic 
forces are then composed of constant value of lift, drag, 
and moment plus a change in these values. These 
forces and moments are resolved into components in 
the coordinate axis directions and are expressed in terms 


of stability derivatives. For example, 


Fy = Xo t+ X,u+ X,uw + Xq (4) 
where 

Xo = force in X direction caused by control move- 
ment 

X, = OX /Ou = rate of change of force in X direc- 
tion due to change in velocity in X direc- 
tion 

X,» = OX /Ow = rate of change of force in X direc- 


tion due to change in velocity in Z direc- 
tion 


In addition to the aerodynamic forces, the change in 
the components of weight in the coordinate directions 
as the attitude of the missile changes must be con- 
sidered. 


Since the mass flow of the jet has been assumed con- 
stant, the jet thrust is constant and does not appear in 
the equations of motion for small perturbations. By 
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using the Durand notation, the equations of motion 


[Eqs. (3)] in the vertical plane become (¢ = 6) 


‘“ X u X w X qd , 
u— u— w+i{— + Wi)¢at+ 
m m m 
Xo 
gcos 0,6 = 
m 
Lu ° Le ; Z, a 
_ u+w— w+(—U, - g+ (5) 
m m m 
3 Zo | 
qymy = 
m m 
M, Me a . a Mo | 
— u— w+g- oa mpv = 
F ; j q+ 79h 7 
where 
m = mM — ct 
[ = mK? 
my = initial mass at time ¢ = 0 (slug) 
ti = mass flow, a constant (slug per sec.) 
K = instantaneous radius of gyration at time / 
(ft.) 
6, = angle of missile axis from horizontal reference 


line for steady-state flight condition 


Jet ForRcES FOR ROCKET-TYPE MISSILES 


In order to determine the jet forces, the summation 
of the effective forces on the particles of the jet must 
be made. Consider the jet shown in Fig. 2. If the 
mass density of the gas of the jet is y and the cross- 
sectional area of the jet is A at any section a distance 
p from the mass center, then for steady continuous 


flow 
yAv = c = constant 
v = ¢/yA 
m, = yA dp 
Therefore 
j lo 
my == of dp = c(ly — |) 
h 


(6) 


j 
me 


II 

in 

a? 
I 

6 

3 


j lo 
>Yomipv = c f pdp = (c/2)(l? — 1,7) 
h 


Hence, in Eqs. (3) the jet 
damping terms are 2Qc(/y — 1;) and Qe(lo? — 1/,”). 
similar manner, the corresponding terms in Eqs. (5) can 
be determined. The equations of motion are, there- 


thrust is cv, and the jet 
Ina 


fore, 
mVx+mVzQ — cw = Fx / 
mV; — mVxQ + 2c(lo —_ 1)Q = Fz (7) 
IQ + c(lo? — h*)Q = Me j 
and the equations in terms of the perturbations 


are 
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X u X w Xq r 
u— u— w+i{—- + Wi}qt+ 
m m m 
Xo 
g cos 6,60 = 
m 
Le ' * ae y Le 
eee eH Oo Pi = q + 
m m m ° 
c ‘ Zo | (8) 
2— (lb —h) q+ g sin 0,6 = 
m m 
M, M,, 14 M, rn 
_ -U— w -- 
? ¥ 7” 7 * 
ane Pe on Mo | 
T (lo “)q= I 


The last set of the equations of motion can be simpli- 
fied by making certain assumptions and approxima- 
tions. 


Thus, for a missile having symmetry about two 
planes and having zero incidence between the longitu- 
dinal axis and the flight path for a steady flight condi- 
tion and flying horizontally at supersonic speed so that 
Z,/m can be neglected in comparison with Uj, the equa- 
tions of motion become 


ua — (X,/m)u + g6 = Xo/m 


: be [ Cc Lo 
w= —w 2 lb —1l)—-—U = 
m L m ( uv | q m (9) 
M,, : 4 F / ° ] 9 M, - Mpo 
I 7 q I (lo 1") 7 ie ] 


Jet Forces FoR AIRJET-TYPE MISSILES 


The jet forces for the airjet-type missile will be de- 
‘veloped in a manner similar to those for the rocket, the 
main difference being that the flow of air from the en- 
trance plane to the fuel injection point must also be 
considered in the case of the airjet. 


The jet terms in Eqs. (3) will be evaluated in two 
parts: one involving the entering air and the other the 
exhaust jet, as shown in Fig. 3. Using the subscript a 
to denote the values pertaining to the airflow and e for 
the exhaust jet, then for steady flow at any cross sec- 
tion 


Va = Ca/ VaAa 
Ce VA. 


I 


where ¢, — Cq = ¢, rate of mass of fuel flow; also 


Mia = YaAqdp 
Mie = YeAdp 


Therefore 


j h mle 
my = af dp + C J dp = Ca(l2 +h) + 
—ls h 


Ce(Lp — ly) 
j "1 vo 
me = caf dv + ¢ [ dv = Cq(v1 — v2) + 
V2 V1 
Ce (Vo — %)f (I 
j 7h “lo | 
> mpu = Ca J pdp + af pdp = 


G i ; C, : : 
(1? eae /,”) + 3 (Io? -_ 1,?) 


where 
v; = velocity of mass flow at point where fuel ; 
added 
v2 = air inlet velocity at entrance plane 
v = jet outlet velocity at exit plane 


The jet thrust can be expressed in terms of the fuel 
mass flow and air mass flow as c(v — v1) + C4(v —% 
Similarly, the jet damping force and moment ar 
respectively, 


2[c(lo os h) ss Calle + 1o)|Q 
[Ca(lo? — h®) + c (lo? — 1210 


The equations of motion for the airjet therefore can 
be expressed as 


mVx + mVz0 — 
[c(vo — V) ca Calte = Vo) = Fy 


mVz — mVxQ + 
2[c(Io <— 1,) +t Ca(lo + 1))|Q = Fy, 


IO + [ca(lo? — le?) + c(lo? — h2)]0 = Me 


or 
. © Aes Ra : 
“ — u— w + (- ‘4, \g + 
m m m / 
Xo 
g cos 0,0 = 
m 
} a : ie ( . “) 
— u+t+w- w+{—U,; — gq + 
m m m ; 
(12 
2 Zo 
[c(Lo = 1;) of Calls oo Io) tie 
m m | 
M, M, i M, 
— u— w - 
I or. se 
l Mo 


] [c(lo? — 1,7) + cally? — 15?) ] g = 


As for the rocket, these last equations for specified 


steady flight conditions can be simplified to 
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i — (X,/m)u + g8 = Xo/m 


eos fu w+ 
Fa m 
oe. 2) + (I, + Jo)] u, | mall 
. , 1 } Ca\l2 0 — -G¢ = 
lm [el ; ‘s* m (13) 
= wt+tgqt 
M)\  — Mo| 


F [e(lo? — hh?) + Ca(l? — 12?)| — rf q I 


SIGNIFICANCE OF JET DAMPING TERMS 


The terms for the jet damping of the rocket given in 
Eys. (9) are 2(c/m )(J) — 11)q for the force in the Z direc- 
tion and (c/J)(/o? — 1:*)q for the pitching moment. It 
is interesting to note that, since /, locates the position 
of the origin of the mass which is eventually ejected 
and Jy locates the pc sition of the detachment of the mass 
from the main body, if the detaching mass originates 
at the point where it is ejected so that /) — /, = 0, then 
there are no jet damping forces. The influence of the 
jet damping force in the Z direction depends on the 
relative magnitude of 2(c/m) (lo — 1,) and U,, which oc- 
cur together in the coefficient for g. Since U, is the 
velocity in the X direction, which in this case has been 
taken as supersonic, it is likely that for many rockets 
U, > (2c/m)(lo — hh), so that (2c/m)(lb — 1) can be 
neglected in comparison with U;. A similar situation 
exists for the terms in the jet damping pitching moment 
where (c/Z)(/)? — 1,°) is compared with \/,/J. In this 
case, however, it is possible to have the value of (c/J) X 
(lo? — 4") 10 to 20 per cent or more of the value of 
M,/I. 

Similarly, for the equations of motion for the airjet, 
the jet damping terms may be small in comparison with 
the values of quantities with which they are combined. 
It should be noted, however, that the magnitude of the 
jet damping terms may be somewhat higher than be- 
fore, since they include the effect of the mass flow of air 
throughout the length of the body. 


SOLUTION OF EQUATIONS 


Even though the magnitude of the jet damping 
terms may be small relative to the other quantities 
with them, the character of the solution of the equa- 
tions of motion for the variable mass missile may be 
quite different from the solution of the equations for the 
constant mass missile because of the fact that in the 
former case the coefficients in the equations are not con- 
stant. 

Since the last two equations of the simplified equa- 
tions of motion for the rocket [Eqs. (9)] and the airjet 
[Eqs. (13)] contain only the two unknowns w and gq, 
it is possible to express them in terms of one unknown 
q. Thus, for the special case of free flight with no 
control forces acting so that X¥y9 = Zp) = Mo = O, both 
sets of equations can be written in the form 


JET- 
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m*j + Amg + (B + Cm)q = 0 (14) 


where, for the rocket, 
A = K~[c(h? — hh? — K*) — M, — 2,K?| 
B= K~*[Z,M, — Zyc(lo? — hh?) + 2Myc(lo — h)] 
C = —M,,U,K- 


and, for the airjet, 
A= K 2[c(lo? —_ 1? — K?) + Ca(lo? = 1,”) md 
M, — ZeK? 
B= K~%\ZyM, — Zw[e(lo? — hi? )+ ca(lo? — l2)?] + 
2M,,[c(lo =e l,) + Calle +- Io)| } 


C = —M,UiK~ 


Changing to m as the independent variable so that 
y = qm, then Eq. (14) may be written as 


m*y" + amy’ + (b + em)y = 0 (15) 


where y’ = dy/dm, etc. The substitution y = mx 
transforms Eq. (15) into Bessel’s differential equation 
in x(m), so that, finally, the solution of Eq. (14) be- 


comes 
g = m~" t+ PIC I, (2Vem) + CY (2Vem)| | 
w= mt 971 C [DJ,(2V/em) + 
E2~V/emJ, « ; (20/em)] + 
C,[DY,(2\/em) + E2/em X 
Y,41 (2 v/em)]} 


(16) 


where, for the rocket, 


» = V(—a + 1)? — 46 
a = —{1+ K-[(h? — 42) — c-(M, 4 
ZK?) } 

b = 1+ K-*{ (lo? — 4?) +e '[—Z,, (lo? —h2)+ 
2M (lo — hi) — (My + ZeK?*)] + 
c~*Z,M,} 

e = —c"K-*M,U;, 

D : M. ‘fe(lo” — 1° + K?) + (cK?/2)(a + 
1—v)-— M) 

E = (cK?/2)M,— 

J = RBessel’s function of the first kind 

Y = Bessel’s function of the second kind 

Ci, Ce = arbitrary constants 


In general, the Bessel’s functions in Eqs. (16) are of 
high order and large argument so that an asymptotic 
approximation can be effectively used. The form used 


herein is 


cos [x — (2/4) — (v7r/2)] 
J x) i - oe 
V (1/2)rx 
e sin [x — (2/4) — (vm/2)] 
Y,(x) =- 7 7 
V (1/2)2x 


The particular problem of the effect of a gust will be 
considered. If we assume that a gust velocity of wo in 
the Z direction is suddenly applied at the time the 





















































202 JOURNAL OF THE AERONAUTICAL SCIENCES~—APRIL, 1950 
0.8 — -——}-- —;—_] 
\ | = 
0.6 —_ | ] T eee 
| | 
04 ‘ N — = --4-—4 ——— 
i XU | | 
0.2 } of}. —}--_J —- 7 AN ! 
we j \ Z N | 
wa 
1 



















































































| 
-0.6 7 
| 
-0.8 —— —— 
“a | | 
ce) 0.08 0.16 0.24 0.32 0.40 048 0.56 064 0.72 0.80 


t (time , sec.) 


Fic. 4. 


missile is in steady free flight, then the boundary con- 
ditions to be satisfied are w = w at f = 0 and gq = Oat 
t = 0. These conditions determine the values of C 
and C, so that Eqs. (16) become 


— [(a/2) 3/ 
w (”.) Ita ~T D a ac 
= : sin 2~/e 
Wo My E2~/ em) 


3/4 
’ ; m , 
(Ym — Vm) + ( ) cos 2V/e X 


Mo 


( V/ my = \/m | 


q 1 m —[(a/2) + (8/4)] 
= ( ) sin 2V/e X 


Wy E2V emy \mo 
/ / | 
(V my — Vm) | 
These values give curves of damped oscillations of 
changing frequency with time. 


COMPARISON OF VARIABLE MASS AND CONSTANT MASS 
ROCKETS 


In order to obtain an idea of the effect of the variation 
in mass on the motion and stability of a rocket, the 
solution for a particular rocket will be compared with 
one for which it is assumed the mass is constant. 

The solution of Eqs. (9), assuming m is a constant and 
c = 0 for the same boundary conditions as before 


namely, w = wy and g = Oat? = 0—is 


w ” a M, ; 
=e | cos wt + _ = | Sin wt 
Wo w mK*w 


= (M,,/mK*w) e™ sin wt 


(18) 
q/ Wo 
where 


a = (1/2mK?)(M, + Z2,K?*) 


Effect of disturbance W 5 on behavior of a variable mass and a constant mass missile. 


l 
—— 


~ ImoK V —K-2 (My + ZyK2)? + 4(ZpMy — My Um) 


Eqs. (17) and (18) have been evaluated for the fol- 
lowing hypothetical values: 


m, = 15 slugs M, = —300 Ibs. ft. sec. 
c = 0.25 slug per sec. K = 2.5 ft. 

Zy = —20Ibs. sec. per ft. h = 5 ft. 

Z, = —20 lbs. sec. i, =0 

M, = —20 Ibs. sec. U, = 1,900 ft. per sec 


The results for the transverse velocity w are shown ii 
Fig. 4. In this figure, the subscript c denotes the con- 
stant mass missile and the subscript v denotes the var: 
able mass missile. It is evident in this case that, since 
the curves have decreasing amplitude with time, both 
the constant mass missile and the variable mass missile 
are dynamically stable. Further, since the amplitude 
of the velocity of the variable mass missile decreases 
more rapidly than that for the constant mass missile, 
the net effect of the jet damping forces and variations 
in mass is to increase the dynamic stability. These 
curves are indicated for too short a time to show the 
effect, but the frequency of the oscillation of the vati- 
able mass missile is increasing with time, whereas the 
frequency for the other missile is constant. 


The ratio of the amplitudes of the pitching velocity 
and transverse velocity, as determined by the envelope 
of the curves of the type in Fig. 4, has been determined, 
and the results are shown in Fig. 5. The amplitude of 
the velocity for the variable mass is always less than for 
the constant mass missile. Since the amplitude of the 


motion damps out rapidly, it is somewhat academic 
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DYNAMIC STABILITY OF 


to compare the ratios of the amplitudes of the motion 
jor a time much over | sec. However, the curves in 
Fig. 5 have been determined for a period of 10 sec. 
to show the trend of the ratios as the time becomes 
large. 

The frequency of oscillation of the variable mass 
missile is not a constant. The ratio of the frequency 
of the variable mass missile and the constant mass 
missile is indicated in Fig. 5. It is evident in this case 
that the frequency of oscillation for the variable mass 
missile only becomes slightly greater than that for the 
constant mass missile even 10 sec. after the disturbance 
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Supersonic Tests of Conventional Control 


Surfaces on a Double-Wedge Airfoil 


MICHAEL PINDZOLA* 
United Aircraft Corporation 


PRESSURE COEF FIGIENT,Go. 


SUMMARY ing of the angle of attack was possible during opera. 

Results of supersonic tunnel tests made to obtain the char tion. 
acteristics of conventional control surfaces are summarized. A A single-pass, double-mirror Schlieren optical system 
symmetrical double-wedge airfoil of 7 per cent thickness, with was used to observe the flow field. Both spark and 


Be $ -r ce ili -e » fi. . ras ~lascs: > Ires- . ° 
15 and 30 per cent trailing edge flaps, was used as a model. Pres time-exposure photographs of the flow pattern were ob. 
sure measurements and Schlieren photographs of the airfoil were 








; : : 3 ; tained. 
taken in a two-dimensional channel at Mach Numbers 1.38, 
1.48, and 1.58. The pressure data were reduced to lift coefficient, 
drag coefficient, lift to drag ratio, section moment coefficient, DESCRIPTION OF MODELS 
flap moment coefficient, and center of pressure points for every 
test angle of attack. Comparisons of the experimental data [he basic model for this investigation was a 6-in, 
were made with the linearized and shock-expansion theories of chord symmetrical double-wedge airfoil of 7 per cent * 
references 1, 2, and 3. Agreement with hones s second- thickness. Detachable flaps were machined to fit the 
ler linear theory and the shock-expansion theory is good except er . me: a : ; 
ee be F cto ieticer airfoil as shown in Fig. 1. Since the sting interrupted 
where such unpredictable phenomena as the following exist: | A ; ‘foi fi 
(1) the shock detachment from the leading edge of the airfoil, the flow over the bottom airfoil surface, flaps were 
(2) the forward movement of the shock emanating at the hinge machined at both positive and negative slopes to obtain 
point of the flap, and (3) separation on the upper flap surface to pressure data for both surfaces. Two sets of flaps inter- 
a limiting case at which the top trailing-edge shock moves to the cepting 15 and 30 per cent of the chord length, respec: 
hinge point. The tests were especially significant in determining ‘ R ; : 
ge} = — Pre ote m8 tively, were made with fixed deflection angles of plus 
the effects of these irregularities on airfoil characteristics. Since ie 50 199° 450 oan ; 
a presentation of all the material would be somewhat lengthy, only @2@ Munuso , 10°, 15°, and 20°. Pressure taps were in- 
representative results are given in this paper. stalled in two rows along the chord with a total of eight 
on the upper forward surface and ten on the flap and \ 
= . portion of the upper rear surface of the wedge. In ad- on 
DESCRIPTION OF TEST EQUIPMENT ae : . : gra 
dition, a check tap was installed on the bottom forward thr 
’ I Sue U.A.C. RESEARCH DEPARTMENT WIND TUNNEL, _ surface. 
in which this investigation was made, is a two- 
dimensional channel of 3.7 by 15.4-in. cross section at 
the test section. Continuous operation is obtained > 
. js > M=1.58 4 
through a vacuum-type system in which part of the > To 6 





mass airflow can be recirculated. This feature allows a sie a | 
simple control of the stagnation air temperature, which | 
was kept constant at 175°F. At this temperature, con- | 
densation was eliminated in the test section. 


Nozzles were designed for Mach Numbers 1.4, 1.5, 30 % DETACHABLE FLAP 





and 1.6 by the method of characteristics, with a bound- 
ary-layer correction of 0.010 in. per in. of the nozzle  M:!5g% 
length downstream of the throat. Subsequent calibra- a 


tion showed the correction to be slightly conservative, 





since actual Mach Numbers of 1.38, 1.48, and 1.58 were , 
obtained. Nozzles were constructed of laminated (b) = 
mahogany with a high luster finish. - 
The model was mounted from the tunnel floor on a bo 
strut of '/4-in. thickness. A biconvex windshield en- thi 
closed this strut and, in addition, the pressure tube leads pl 
and a rod to set the airfoil angle of attack. The chang- co 
Presented at the Aerodynamics Session, Annual Summer Meet- — 
ing, I.A.S., Los Angeles, July 21—22, 1949. ch 
an 
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METHOD OF OBSERVATION 


A complete set of pressure data and a Schlieren photo- 
graph were taken at angle-of-attack increments of 1° 
through a range indicated in Table 1. 


TABLE 1 


— Flap Angle-of-Attack Range 

% Chord Deflection, Deg. M =1.838 M=1.48 M =1.58 
15 0 —3to3 -—T7to7 —9to9 
15 5 —3tol -—7to4 —9tob 
15 10 —3tol —7to3 —Sto 5 
15 15 —3tol -—7to3 —S to 5 
15 20 —4to0 —7to2 —9to4 
30 0 —3to3 —7to7 —9to9 
30 5 —4to3 -—7to4 —9tob 
30) 10 —4tol -—7to3 —9to4 
30 15 —7 tol —9to2 
30 20 —T7Tto-1 —-9tol 


Limits were dictated by bow shock detachment at the 
negative angles and by either complete separation or 
hinge shock detachment at the positive angles. 

The interruption of flow caused by the sting on the 
bottom surface of the airfoil made pressure readings on 
this surface meaningless. Therefore, to obtain a com- 
plete airfoil pressure distribution for any given set of 
conditions, data from two separate runs were used. For 
example, Fig. la shows the airfoil with a 30 per cent 
chord, 10° downward deflected flap operating at a = 6 
1.58. 


and .\/ The pressure distribution data for 


Separation effects at 1/J = 1.48 on 15 per cent chord flap deflected 15°. 


surface A are obtained from the run on the configura- 
tion in Fig. lb at a = 6°. Corresponding data for sur- 
face B are obtained from the run on the configuration in 


Fig. le ata = —6°. 


DISCUSSION 


Since the nature and effects of the irregularities in the 
flow pattern cannot be predicted by existing theory, a 
discussion of them is of interest. This is most easily 
accomplished with the aid of Schlieren photographs and 
corresponding pressure distribution diagrams shown in 
Figs. 2, 3, and 4. These typical examples were taken 
from the runs at Mach Number 1.48 on the 15 per 
cent chord flaps. The theoretical results shown in 
the diagrams were computed by the shock-expansion 
theory. 


Separation 

The series of diagrams and photographs of Fig. 2 
shows various stages of boundary-layer separation. As 
separation increases, the pressure on the portion of the 
flap affected approaches free-stream static pressure. 
Where separation is complete, as at a = 3°, the trailing- 
Lift 
and pressure drag of the airfoil are decreased with the 
The center of pressure of 


edge shock moves forward to the hinge point. 


existence of this condition. 
the airfoil moves forward. 
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Hinge Shock Forward 

The plots and photographs depicting the hinge shock 
forward irregularity are shown in Fig. 3. Note that the 
flaps in the photographs are deflected upward, simulat- 
ing the bottom surface of an airfoil with downward flap 
deflection. Under normal conditions, a shock originates 
at the lower surface hinge point of the airfoil as shown 
in Fig. 3 at a = —2°. This shock turns the air stream 
through the angle of the flap deflection. However, 
when previous disturbances slow the air stream down 
to a Mach Number at which this shock detaches, it 
moves forward, in effect making a longer flap. This 
condition is shown at a = —4°. The increase in pres- 
sure behind this shock thus acts over a greater portion 
of the airfoil. Although it would appear that the ex- 
perimental pressure coefficients on this effective flap 
would drop below the theoretical, such was not the case 
on any flap where the hinge shock moved forward. 
An increase in lift and pressure drag results, in addition 
to a rearward shift of the location of the center of pres- 


sure. 


Detached Bow Shock 


Comparative plots and photographs at points of 
attached and detached bow shocks are shown in Fig. 4. 
At a = —7°, the upper forward surface is inclined at 
too great an angle for the particular free-stream Mach 
Number, and the bow wave separates from the wedge 
tip. This causes a high pressure at the leading edge of 
the airfoil which decreases steadily to the theoretically 
predicted value. The results are an increase in lift 
curve slope, an increase in drag, and a forward shift of 
the airfoil center of pressure location. 


COMPARISON WITH THEORY 


The experimental results of this test were compared 
with the first- and second-order linear theories of ref- 
erences | and 2, as well as the shock-expansion theory 
of reference 3. In all cases, the results agree favorably 
with the second-order linear and _ shock-expansion 
theories up to points where irregularities occur. A rep 
resentative set of comparative characteristic curves 
is given in Figs. 5 through 9. These were obtained at a 
Mach Number of 1.58 on the airfoil with a 30 per cent 
chord deflected 15°. During test, the bow shock was 
detached at angles of attack of —S° and —9° for this 
run. Partial separation of the flow on the top flap 
surface occurred at 0° angle of attack. Separation was 
complete at the 1° and 2° points. 


Lirt COEFFICIENT 


The curves of the lift coefficient, C,, versus angle of 
attack, a, are shown in Fig. 5. The experimental curve 
crosses the zero ordinate at —3.5°. Corresponding 
values for the first-order linear, second-order linear, and 


shock-expansion theories are —4.5°, —3.7°, and —3.9°, 
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Fic. 7. Experiment vs. theory at J = 1.58 on airfoil with 30 
per cent chord flap deflected 15°. 


respectively. In the region of zero lift, the experimental 
curve of lift coefficient is practically linear. However, 
when separation occurs, a loss in lift results. In passing 
from an attached to a detached bow shock condition, the 
slope of the lift coefficient curve increases rapidly. 
However, once the shock is detached, the coefficient 


increases slowly. 


Drag Coefficient 

Fig. 6 shows the drag coefficient, Cp, curves obtained 
by experiment and theory. The experimental curve 
follows the familiar parabolic shape for angles at which 
the lift coefficient curve is linear. Agreement with the 
second-order theory is good in this region. However, 
when irregularities occur, deviations from the parabola 
exist. In the case of separation, the drag coefficient de- 
creases, while with a detached bow shock the op- 
posite is true. These drag coefficient values represent 
pressure drag only, and friction drag may alter the 
shape of the curve, especially in the region of sepa- 


ration. 
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Center of Pressure 


The center-of-pressure results are shown in Fig, 7 in 
terms of the per cent chord from the leading edge 
Agreement of the second-order linear and_ shock-ex 
pansion theories with experiment is again good. At the 
angle of zero lift, the center of pressure is at infinity 
As the angle of attack is changed in either direction 
from the zero lift angle, the center of pressure ap- 
proaches the mid-chord position. 


Section Moment Coefficient 


Fig. 8 shows the section moment coeflicient, Cy 
curves. These moment coefficients were obtained by 
integrating the elemental moments about the 50 per 
cent chord position. A positive moment was taken as 
that tending to raise the leading edge of the airfoil 
Agreement of the experimental results is best with the 
second-order linear theory. The effect of irregularities 
in the flow pattern was to increase the slope of the mo- 


ment coefficient curve. 


Flap Moment Coefficient 


The curves of flap moment coefficient, Cy, are shown 
in Fig. 9. These coefficients are for the flap section 
about the hinge point and are useful in determining 
control loading. The second-order linear theory ap- 
proximates these loads accurately up to a point where 
separation on the flap occurs. Since separation tends to 
increase the pressure on the upper flap surface, the flap 
hinge moments decrease. 
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SUMMARY CURVES 6 7 T a 


Figs. 10 through 14 show the effect of flap deflection 
angle, 6, on the displacement of the characteristic i 
curves. The plots were made from experimental data of 
runs at Mach Number 1.58 on the 30 per cent chord 
flaps which include the run compared with theory 2 
above. The displacement of the curves is consistent in 
regions where no irregularities exist. However, where 











irregular conditions are present, the curves deviate in 0 
accordance with the irregularity as discussed previously. <. 
The deviations at angles of attack of —8° to —9° are &@ 
caused by a detached bow shock in allruns. Partial sep- = = “ 
aration has little effect on all but the flap moment coeffi- 8 vs 4 
cients. Complete separation has an appreciable effect 5 specil 
on all characteristics as shown by the plots of the 15° and -.4 -+ cat. the t 
20° deflection runs at the high angles of attack. These ial casita how 
effects are canceled somewhat on the 20° run by the wc : (SHOCK EXPANSION) = 
existence of the hinge shock forward condition. Similar val Me Senne t —————; ' seem ope 
. . . ~ | 2 HINGE SHOCK FORWARD . 
trends are obtained with either 15 or 30 per cent chord | . ae, cei 
flaps at other Mach Numbers, although the displace- 4 COMPLETE SEPARATION 
ments of the curves are less with the shorter flap. -.8 ~ : i = Ke 
Fig. 15, a cross-plot of Fig. 10, shows the effect of ANGLE OF FLAP DEFLECTION, 3 - OEGREES T 
flap deflection angle, 6, on the lift coefficient at constant fFy¢. 15. Effect of flap deflection angle, 6, on Cr at constant 
angles of attack. Theoretical curves computed by the angle of attack, a, and M = 1.58 with 30 per cent chord flap, : , 
snee 


shock-expansion method are included to indicate the 

i esis ; ; struc 
effectiveness of the flaps. At any constant angle of sults indicate less increase than the theoretical, since the TI 
attack, the lift coefficient increases with an increase in experimental pressure distributions on the flap are 





P ‘ - ° may 
flap deflection angle. However, the experimental re- slightly lower than those predicted. ; 
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Minimum Weight of Stiffened Cylindrical 
Shells in Pure Bending 


W. R. MICKS* 
The RAND Corporation 


SUMMARY 


Methods previously developed for determining minimum 
weight of stiffened cylindrical shells in bending! are extended 
to include the weight of the circumferential frames. For a 
specified diameter and bending moment it is possible to minimize 
the total weight of shell and frames. Curves are presented 
showing the variation of frame weight, panel weight, and total 
weight with variations in structural index and frame spacing. 
These curves may be used to indicate the weight penalty incurred 


by using nonoptimum design. 


DERIVATION 


_— TOTAL WEIGHT of the stiffened shell is con- 
sidered to be composed of two parts: the weight 
of the circumferential frames and the weight of the 
sheet-stringer combination that comprises the panel 
structure between frames. 

The total average weight per inch length of shell 


may be written 


W,’ = W,’ + Wy’ (1) 
where 
Wp’ = average panel weight per inch length (prime 
indicates ‘‘per inch’’) 
W,’ = average frame weight per inch length 


Assuming that the cross section is circular and that 
the panel proportions are the same all around the 
circumference, the panel weight is given by 


Wp’ = rDiw, (2) 
where 
D = diameter of shell 
/ = average thickness of panel 
w, = density of panel material 
The weight of frames (per inch length of shell) is 
Wr’ = xDAyur/L (3) 
where 
A, = cross-sectional area of the frame 
L = frame spacing 
we = density of frame material 


It is shown in reference 2 that the stiffness required 
of the frame cross section, in order to prevent general 


instability, can be expressed by 
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Kpl p = C,MD? i (4) 
where 
Ey = modulus of elasticity of frame material 
I], = moment of inertia of the frame cross section 
C, = a dimensionless coefficient? 


An analysis of test data has shown that a stiffened 
having a value of C;, above a critical value 
(actually a critical range) will exhibit a panel type of 
failure wherein the buckles between frames, 
leaving the frames undistorted. As the present analysis 
is based on the assumption of designing to prevent a 
general instability failure, the value of C; will be chosen 


cylinder 


shell 


accordingly. 

Assuming that the value of / is known, the required 
moment of inertia, J,, must be related to the weight 
(or cross-sectional area) of the frame. This will be 
accomplished by use of the dimensionless shape param- 
eter! ky, where ky = /,/A,’. 

This parameter may be thought of as an “efficiency 
factor’’ for the frame cross section in bending. 

Transposing, 


Ip = biAy? (5) 


Substituting Eq. (5) in Eq. (4), the required cross- 
sectional area of the frame is 


A,*? = C,MD?*/kiLE, (6) 


This may be written 


(DM""/L'"")(C,/ksEpr)”* (Ga) 


Ay = 
Substituting Eq. (6a) in Eq. (3), the average frame 
weight per inch becomes 


Wy’ 


= (4D? M*wy/L’)(C,/kiEp) (7) 

It can be seen from Eq. (7) that, as C, decreases’ 
the frame weight will decrease. However, a lower 
limit is fixed for C, by the stipulation of a panel-type 
failure. For any value of C, above this limit, the 
frames will obviously be overweight, as far as general 
instability is concerned. 

It should be noted that Eq. (7) shows the frame 
weight to be proportional to the factor wp/WEp, 
which is density divided by the square root of the 
modulus of elasticity. This indicates that, with other 
factors remaining the same, the use of some material 
other than aluminum alloy may effect a saving in frame 
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weight. For example, the factor wy/+7/ Hp is 2.55 X 
10~-° for magnesium alloy as compared with 3.09 X 
10~° for aluminum alloy. Therefore, the ratio of the 
weight of a magnesium-alloy frame to the weight of an 
aluminum-alloy frame having the same cross-section 
proportions and the same value of // is approximately 
0.826. It should be that, 
although the above figures show that a lighter frame 


pointed out, however, 
could be made from magnesium, the magnesium frame 
will be larger in cross section. In order to have the 
same EJ as the aluminum frame, all dimensions of the 
magnesium frame cross section will be multiplied by the 
factor [E(Al)/E(Mg)]"‘, which equals 1.13. There- 
fore, the amount of weight that can be saved is de- 
pendent on any size restrictions that might be imposed 


4 


on frame depth or other frame dimensions. 
Eqs. (7) and (2) are added to give the total average 
weight per inch length. 


P , rD?\' ( C; ) , .e 
W,’ = rDiu, (8) 
es” ae , 


Assuming that the frames and panel are made from 
materials of the same density, w, = wp. 
Also let 


C2 = (C;/RksEv)” (9) 


Eq. (8) may be expressed in terms of total average 
cross-sectional area. (NOTE: The term “total average 
area”’ will be used to denote average volume of panel 
and frames per unit length of cylinder.) 


(10) 


A = rDi + (rQD?M/L") 


The average panel thickness / may be expressed as 


E = Gmaz./fe = 4M/eD7, (11) 

where 
maz. = Maximum compressive load per inch width 
te = allowable compressive stress (average) in 


the panel 
Substituting Eq. (11) in Eq. (10) 


A = (4M/Df.) + (#C,.D°?M"*/L””) (12) 


In order to compare different configurations (different 
ratios of L/D or different types of panel construction) 
of stiffened cylinders at different intensities of loading, 
Eq. (12) will be expressed in terms of the structural 
index .//D* for members in bending and the non- 
dimensional frame spacing L/D. 

By algebraic manipulation, Eq. (12) is put in the 
form 


A — 4(M/D*)’ Cs 


: so (13) 
M*” iP r (if /D*)' (L/D) 


In this form, the equation eliminates all effects of 
size and considers only proportions and stresses. 
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When any two geometrically similar cylinders com. 
posed of the same material are loaded in bending ty 
equal values of 1//D*, they will develop equal valyes 
of A/M”’, regardless of their relative sizes. If \ 
represents the ultimate moment, then A is the tot) 
average area of cross section required to carry this 
A/M™ is a 


measure of the weight of cylinder required to carry , 


value of moment. Therefore, direct 
given value of moment at a specified value of J//D 
(NOTE: 
caused by 
changes of the stress-strain curve with variation ; 


This neglects possible second-order effects 
change of stress gradient with size and 


size.) 

Another advantage of expressing Eq. (13) in terms oj 
the structural index is that 1//D*, together with L/D. 
can be converted into the structural index g/L for 
loaded 
available curves of allowable average stress f, ys, 


axially compression panels. As_ there are 
q/In for many types of sheet-stiffener panels,' data 
from these curves may be used directly in estimating 
the weight of a stiffened cylindrical shell. The value 
of g/Lo is determined as follows: 

The load per inch width g is taken as the maximum 


value for the cylinder in bending. 


Q = mar. = 4+M/xD? 


Assuming end restraint to be zero, the effective panel 


length Ly is equal to the frame spacing L. Therefore, 


qd Lo = 4\/ rD?L (14 
Expressed in terms of the structural index, 
g/Lo = 4(M/D*)/2r(L/D) 


After g/L» is calculated, the value for f, in Eq. (13 
is read from a curve of f, vs. g/Lo corresponding to the 
type of panel used in construction of the shell. Such 
a curve is shown in Fig. 5. 


EXAMPLE 


To illustrate the methods outlined above, Eq. (1s 
will be plotted for an example case, as shown in Figs. 
1-3. 

The following assumptions are made: 

(1) The lower limit of C; is 1/16,000.° 

(2) The value of the coefficient C: remains constant. 
This implies that, after a frame material is specified 
and a value is chosen for C,, the cross sections of all 
frames are geometrically similar (or, more generally, 
have equal values of 4). 

(3) The panel is constructed of 24ST sheet with 
formed stiffeners designed to optimum proportions 
Values for the allowable compressive stress f, will be 
read from the curve of f, vs. g/Lo shown in Fig. 9 
This curve, reproduced from reference 1, is an envelopé 
curve representing the upper limit of stress which has 
been attained in tests of 24ST sheet-stiffener panels 


having formed stiffeners. 
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The value of the shape factor ky is based on the pro- 
portions of the frame cross section shown below. 
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SHAPE OF FRAME CROSS SECTION 


Assuming Hy = 10° lbs. per sq.in. and substituting 
the above value of ky in Eq. (9), ¢ 
109 X 10-*, 


In evaluating . 


» is found to be 


A/M 


is held constant as \//D® is varied. 


from Eq. (13), the ratio L/D 
A value of g/Lo 
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is calculated from Eq. (14a) for each combination of 
L/D and M/D*, and the corresponding value of f, is 
read from Fig. 5. Then A/M”* is evaluated, and the 
point is plotted on the curve of A/M™” vs. M/D*. 
This procedure is carried out for the values of L/D 
shown in Fig. 1, the range of values being chosen so 
as to include the optimum value of frame spacing. 
The curves of Fig. 1 show directly the variation in 
total average area for different values of structural 
index and frame spacing. The cube root of the struc- 
tural index was used for purposes of contracting the 
horizontal scale. 

It is probable that in the design of a fuselage both the 
diameter and ultimate moment will be specified, thus 
fixing the value of \//D*. It remains then to deter- 
mine the optimum value of L/D for the specified value 
of \//D*. For the case represented in Fig. 1, there 
appears to be one optimum value of L/D for the entire 
range of 1//D*—namely, a value of about 0.10. 

To determine more accurately the optimum value of 
L/D, the curves of Fig. 1 are cross-plotted as shown 
in Fig. 2. The optimum value of L/D is seen to be 





AND PANEL PLOTTED SEPARATELY 


approximately 0.115 for the range of 1J/D* shown. As 
the value of L/D is decreased from 0.4 to the optimum 
value, the curves of Fig. 2 show a gradual decrease in 
A/M**. However, any further decrease in L/D 
results in a sharp rise in A/M”’, indicating that, if 
nonoptimum frame spacing is used, the greatest weight 
penalty is likely to be on the side of closer-than- 
optimum frame spacing. 


The range of panel compressive stress f, for the 
curves of Fig. 1 is from 10,000 (left end) to 50,000 Ibs. 
per sq.in. (right end). It is interesting to note that 
the stress corresponding to the optimum 1//D* for 
nearly all values of L/D lies in a range of from 26,000 
to 32,000 Ibs. per sq.in., indicating that the stress for 
optimum design is considerably lower than the maxi- 
mum allowable stress for the material used. 


Assuming that 50,000 Ibs. per sq.in. is the maximum 
allowable stress for the 24ST panel, the dot-dash 
portion of the curve for L/D = 0.10 (Fig. 1) shows the 
effect of holding the stress constant at 50,000 Ibs. per 


sq.in. while increasing 1//D*. 
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The first and second terms of Eq. (13) represent, 
respectively, the amount of panel material and the 
amount of frame material. Although it is the sum of 
these terms which must be minimized to give minimum 
weight, it is desirable to know the relative amount of 
weight in the panel material and in the frames. Fig. 3 
shows the values of A/M" plotted separately for the 
panels and frames. For the near-optimum value of 
L/D = 0.10, it can be seen that the amount of material 
in the panel is roughly from three to six times that in 
the frames, depending on the value of //D*. How- 
ever, as the frames are spaced closer than optimum, the 
ratio of frame material to panel material rises sharply. 

Fig. 4 shows the result of using a less efficient frame 
cross section, the value of /,/A,* being one-tenth that 
of the first example. As might be expected, the use of 
heavier (for a given value of FJ) frames causes the 
optimum spacing to increase over that indicated for 
the first case, and the total area required is likewise 
increased. As the same panel construction was used, 
the values of the first term of Eq. (13) remain the same 
as for the first example. It was necessary to calcu- 
late new values only for the second term of Eq. 
(13) in order to plot the new set of curves of A, M*“* vs. 
M/D*. 
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CONCLUSIONS 


As this analysis considers only stresses and propor. 
tions, it does not include possible effects of size limita. 
tions such as minimum gages, minimum size of stringers 
and frames, etc. These limitations must, of cours. 
be considered in any actual design, since it is likeh 
that optimum proportions could not be maintaine; 
throughout the structure. 

The value of C, chosen for the example is for a circylay 
cross section. 
other shapes. 

It should be noted that the frames in the aboy 


A different value would be needed fo; 


analysis are designed to prevent general instability 
when the shell is subjected to pure bending. Other 
design -such as transverse shear, local 
loadings, and pressure loads—may be critical for the 


conditions 
frames in an actual design. 
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Analysis and Design of Stiffened Shear 
Webs 


PAUL H. DENKE* 
Douglas Aircraft Company, Inc. 


ABSTRACT 


A simplified graphical solution of a previously published elastic 
energy analysis of shear panels is presented. Comparison is 
made between the results of the analysis and experimental data, 
and satisfactory agreement is found. A qualitative description 
of the failure of flat and curved shear panels is given, together 
with an example of the stress analysis of flat panels. It is con- 
cluded that a rational, easily applied, and satisfactorily accurate 
method has been developed for determining the stress distribu 
tion in stiffened shear webs. It is expected that future develop- 
ment, some of which has already taken place, will permit the pre- 
diction of ultimate loads not only for thin web beams but also for 


curved panels typical of wing and fuselage structure. 


INTRODUCTION 


hee PRESENT PAPER provides a simplified graphical 
solution of a previously published shear panel 
analysis.' It also gives a comparison with test data and 
shows some of the applications in stress analysis and de- 
sign. The previous paper presented a mathematical 
solution of the problem, but the results were not in a 
form suitable for practical calculation. This difficulty 
has been remedied through the discovery of a means of 
graphical presentation of the important variables. 


The experimental comparison is with the N.A.C.A. 
data on diagonal tension beams.? Other experimental 
studies of the stress distribution in shear panels are 
available, but time did not permit a critical comparison 
with all. Agreement with test data is satisfactory in 
most cases, and such discrepancies as exist can probably 
be explained on the basis of difficulties involved in the 


testing technique. 


The section on engineering applications is necessarily 
limited for lack of space, but it is hoped that the discus- 
sion given will serve as the basis for further develop- 
ment. It is deemed advisable, at present, to avoid par- 
ticularizing too much with regard to the estimation of 
ultimate shear loads, because these depend on the detail 
design of the panel which varies a great deal among the 
various air-frame manufacturers. Greatest emphasis is 
therefore laid on the determination of the stresses exist- 
ing in various portions of the panel for a given shear 
load. Whether or not these stresses will cause failure 
in any of the panel components requires additional cor- 
relation between analysis and test results. 

Presented at the Structures Session, Annual Summer Meeting, 
1.A.S., Los Angeles, July 21-22, 1949. 

* Structures Engineer. 


Historical Survey 


The following brief outline is not meant to be exhaus- 
tive but is intended only to trace the growth of the prin- 
cipal methods of shear panel analysis so far proposed. 

Immediately after the publication of Wagner’s class- 
ical treatise on the pure tension field,* it became clear 
that the assumption of zero diagonal compressive stress 
was entirely too conservative, even though it did pro- 
vide a clear qualitative understanding of the behavior of 
thin-web beams in shear. In order to avoid this dif- 
ficulty, three principal methods of analysis or design 
arose, of which the first was the purely experimental ap- 
proach whereby components such as spar webs, wing 
panels, or fuselage panels were tested to destruction. 
This is the soundest means of determining the strength 
of a shear panel, provided suitable testing methods are 
employed. However, it is expensive and time-consum- 
ing and is not well adapted to certain purposes, such as 
the systematic search for optimum panel proportions, 
because of the large number of variables involved. 
Moreover, the experimental inaccuracies that inevitably 
creep in make it advisable to have on hand an analysis 
method that can serve as the basis for interpreting test 
results. 

The second approach was through the modification of 
the Wagner equations. This was accomplished most 
simply by assuming that the principal diagonal com 
pressive stress was not zero but constant and equal to 
its value at the critical buckling load. 
type have been developed by Langhaar,* Sechler and 
Dunn,*® Niles and Newell,® Bruhn,’ and others. 
methods also tend to be conservative because many 


Analyses of this 
These 


webs buckle at a negligible fraction of their ultimate 
strengths, and the computed diagonal compressive 
stresses, assumed constant after buckling, are therefore 
too small. 

A more realistic modification of the Wagner equations 
was proposed by Kuhn.* In his analysis, the fraction of 
the shear load carried in diagonal tension is assumed to 
be given by an empirically determined function of the 
The remainder of the load is carried as in 
Since the various empirical con 


shear stress. 
a shear-resistant web. 
stants have been evaluated to fit experimental results, 
the analysis can be expected to be reliable in the range 
covered by the test data, provided suitable testing tech 
niques have been employed. The method is not well 
adapted to the analysis of heavy webs, since it provides 
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Fic. 1. Average web longitudinal tensile stress. 


no means of computing the bending loads imposed on 
the stiffeners by the buckled sheet. 

The third approach to the shear panel problem was 
simply to make use of established mathematical pro- 
cedures to determine a reasonably close approximation 
to the true distribution of stress existing in the web and 
its stiffening elements. The mathematical instrument 
usually employed was the principle of virtual displace- 
ments, or minimum potential energy, various forms of 
which were developed by Rayleigh, Ritz, and Galerkin. 
These methods make use of the well-known expression 
for the strain energy stored in the buckled web.'° The 
first analysis of this type appears to have been that of 
Kromm and Marguerre,’ who applied the energy prin- 
ciple to the problem of the infinitely long shear panel 
subjected to stresses greater than critical, for the two 
cases of infinite and negligible transverse and longitudi- 
nal stiffening. The mathematical basis of the present 
paper, first presented in 1943, derives from the work of 
Kromm and Marguerre, although certain points of pro- 
cedure are quite different. Recently, Levy and others 
have published analyses employing the energy method 
or the equivalent differential equations of von Karman. 

Besides the shear panel analysis, the writer has on 
hand certain unpublished investigations into the be- 
havior of long, flat, and curved panels subjected to com- 
bined shear and compression. 


Nomenclature 


In the following discussion, “‘longitudinal’’ and 


“transverse”’ refer, respectively, to the long and short 


dimensions of the panel. Thus, web stiffeners on a 
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thin-web beam are longitudinal stiffeners if they ar 
closer together than the flanges; otherwise, they ar 
transverse. On a fuselage panel the longerons are orqj. 


narily closer together than the frames. 


The longerons 


are therefore longitudinal stiffeners, and the frames are 
transverse stiffeners. 


NOTATION 


The notation of reference 1 has been revised as follows: 


W 
L 


Wr 


Ri 


ll 


ll 


ll 


ll 


ll 


Il 


Fic. 





panel width free to buckle, in. 

panel length free to buckle, in. 

web thickness, in. 

average web shear stress, Ibs. per sq.in. 

effective web modulus of elasticity, lbs. per sq.in. 

effective longitudinal stiffener modulus of elasticity, 
Ibs. per sq.in. 

effective transverse stiffener modulus of elasticity 
Ibs. per sq.in. 

longitudinal stiffener cross-sectional area, sq.in 

transverse stiffener cross-sectional area, sq.in. 

longitudinal stiffener moment of inertia, in.‘ 

transverse stiffener moment of inertia, in.‘ 

longitudinal stiffening ratio 

transverse stiffening ratio 

web average longitudinal tensile stress, lbs. per sq.in. 

web average transverse tensile stress, lbs. per sq.in 

web average diagonal tensile stress, Ibs. per sq.in. 


web average diagonal compressive stress, Ibs. per sq.in, 


longitudinal stiffener local bending load, lbs. per in. 
transverse stiffener local bending load, Ibs. per in. 


longitudinal stiffener local bending moment, in.lbs. 


transverse stiffener local bending moment, in.Ibs. 
panel secant shear modulus = f,/y,, lbs. per sq.in. 


average panel shear strain 


lO ; 16 





Ke, 


2. Longitudinal tensile stress correction factor 
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DESIGN OF 


mam 
we 





Fic. 3. 


Average web transverse tensile stress. 


Ar, = effective longitudinal stiffener area, sq.in. 

p = stiffener radius of gyration, in. 

GRAPHICAL PRESENTATION OF THE SHEAR 
PANEL ANALYSIS 


The graphical method discussed in this section was 
developed in collaboration with D. M. Conner, of Doug- 
las Aircraft. Mr. Connor also contributed to the sec- 
tion on the qualitative behavior of shear panels. Ac- 
knowledgment is likewise made to S. Noskoff from 
whose computations the design charts were constructed. 

The four basic independent dimensionless groups in- 
volved in the analysis are the stress ratio f,WL/(E?’*), 
the panel aspect ratio L/W, and the transverse and 
longitudinal stiffening ratios Ry and R,. The effect of 
the quantity Rg of reference 1, which relates to the 
bending flexibility of beam flanges, is accounted for in a 
later paragraph by modifying Rr and R, in order to 
avoid introducing a fifth or sixth independent variable. 
It is therefore necessary to plot the important panel 
stresses, such as the diagonal tensile stress 6, as func- 
tions of four independent variables. In general, these 
are too many to plot, and the mathematics gives no clue 
as to how matters can be simplified. 

However, it is found by calculation that the various 
panel stresses can be plotted as functions of the stress 
ratio and the panel aspect ratio for Rr = 1 and R; = 
l. The effect of the actual panel stiffening can then be 
accounted for by means of a multiplicative factor K, 
which is a function of Ry and R,. In all cases, the ratio 
of the desired stress to some power of the shear stress is 
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plotted rather than the stress itself, in order that the 
curves shall have asymptotes (see Figs. 1 to 17). Asa 
result, it is possible to estimate values closely, even when 
they run off the curve sheet, and this happens only 
rarely because of the logarithmic scale for f,VL/(£t??*). 
The associated K factor is simply the ratio of the asymp- 
totic value of the function for the given stiffening ratios 
R,and R, to its asymptotic value for Rr = land R, = 
1. An example of the use of these charts is given in a 
later section. 

The following quantities are presented graphically: 
fu, fr, 1, Wr, Wr, Mz, Mr, and G,. Most of the im- 
portant panel loads and stresses are thus represented, 
and the remaining ones can be obtained by straightfor- 
ward calculation. It was not found possible to plot the 
depth, length, and angle of the web buckles in simple 
fashion, but these are not ordinarily required in analy- 
sis. However, further investigation may disclose a 
means of plotting these quantities also. 

The charts are based on a value of Poisson’s ratio 
equal to 0.33. 

It has been found that, by comparison with calcula- 
tions made from the original analysis, the panel stresses 
obtained graphically are given to slide-rule accuracy. 
The charts therefore comprise a solution of the mathe- 
matical analysis complete and accurate enough for en- 
gineering purposes. 


COMPARISON WITH TEST DATA 


The stiffener stresses for the beams of series I, II, and 
III of reference 2 were computed and compared with 
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the test results, as shown in Figs. 18, 19, and 20. The 
method of computation is similar to the example givey 
in a later paragraph. However, in the case of beams 
having stiffeners on one side of the sheet only (these are 
denoted by ‘‘S’” in the beam code number), it was nec. 
essary to account for the fact that test stresses were 
measured only on the leg of the stiffener attached to the 
sheet, so that bending stresses were not eliminated. 
The correction was made according to the simple theory 
of flexure. 

It should be noted that the computed stiffener stress 
is due only to the axial load introduced by the buckled 
web plus the bending due to the fact that the web is ec. 
centric with respect to the stiffener. However, the 
range of test scatter indicates that extraneous stresses 
have affected the measurements. These additional 
stresses may be due to stiffener column action, stiffener 
torsional instability, or lateral loads imposed on the 
stiffener by the buckled web. However, the method of 
the present paper indicates that the web wrinkles pro- 
duce small bending moments in the stiffener. If it be 
assumed that these moments are resisted only by the leg 
of the stiffener attached to the web, it is still only pos. 
sible to account for a portion of the observed scatter. 

In an effort to eliminate the extraneous stresses, the 
strain gage readings were averaged up and down the 
length of the stiffener.?. But it does not follow that the 
resulting average would represent the sum of the com- 
pressive plus eccentric bending stress, as supposed. 
Thus, it is not surprising that the largest discrepancies 
between computed and test results occur for beams with 
single web stiffeners. (Beams having stiffeners on both 
sides of the sheet are denoted by ‘“D”’ in the code num- 
ber.) 

It should be kept in mind in examining Figs. 18, 19, 
and 20 that no empirically determined functions have 
entered the analysis at this point, except the experi- 
mentally determined variation of effective web modulus 
with shear stress discussed later. 


ANALYSIS AND DESIGN 


Failure of Flat and Curved Shear Panels 


This section is a qualitative description of the be- 
havior of both flat and curved shear panels at failure. 
It is believed that the present analysis may be applied, 
with suitable modifications, to curved panels in any 
case where the stiffening elements are sufficiently close 
together so that the sheet bounded by them is, for all 
practical purposes, flat. For example, a_ fuselage 
panel having a radius of curvature of 60 in., with lon- 
gerons 6 in. on center, should come under this classi- 
fication. 


Web Failure 


Shear panel webs normally fail through diagonal 
tension. Although basically simple, the mechanism 0! 
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diagonal tension is complicated by the following factors: 
(1) The diagonal tensile stress is not uniform but, in 
the ideal web, tends to be maximum at the center and 
minimum at the edges. (2) The effect of wrinkling is 
to introduce bending stresses that add to the diagonal 
tension. These stresses likewise vary throughout the 
web and, in the ideal case, tend to a maximum at the 
center. (3) In the actual case, the riveted attachment 
of web to stiffeners produces stress concentrations that 
cause failure to start at the panel edges. Fig. 21 
shows a diagonal tension failure. 
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Fic. 16. Panel secant shear modulus. 
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Longitudinal and transverse local bending load and moment correction factors. 


Stiffener Failure 

The term “‘stiffener’’ is here used in a general sense 
and may denote the stiffeners or flanges of a tension 
field beam, the stringers or ribs of a wing panel, or the 
longerons or frames of a fuselage panel. The following 
types of loads contribute to stiffener failure : 

(1) Axial Loads.—In the case of a beam, the longi- 
tudinal tensile stress f, is transmitted into the stiffener, 
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Panel secant shear modulus correction factor. 
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Fic. 19. Stiffener stress vs. shear load. N.A.C.A. beams series II. 
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producing the compression load P;,, as shown in Fig. 22. 
A stiffener that failed locally in axial compression is 
The transverse tensile stress fy pro- 
In the case of a 


shown in Fig. 23. 
duces similar loads P; in the flanges. 
fuselage panel, the transverse tensile stress fy has an 
inward component because of curvature, as a result of 
which radial loads per inch IV; are produced on the 
stringers (see Figs. 24 and 25), which in turn introduce 
compressive loads P, in the frame. 

(2) Eccentric Bending Load.— Besides producing com- 
pressive stress in the stiffener, the load P; mentioned in 
reference | also introduces bending if the centroid of the 
stiffener does not lie on the skin line. Thus, in Fig. 26, 
the stiffener bending moment is P72. 

(3) Column Action and Torsional Instability. 
stiffeners may have a tendency to fail as columns under 


Beam 





Stiffener stress vs. shear load. 


N.A.C.A. beams series III. 

the action of the stiffener axial loads. However, there 
are certain characteristics of shear panel behavior 
tending to minimize the likelihood of column failure, as 
is demonstrated by the following argument, which is 
thought to be valid in the case where the stiffeners are 
much closer together than the flanges. 

Consider the curved stiffener shown in Fig. 27. The 
curvature may either be built in or be due to the bend- 
ing effect of the various bending loads, discussed else- 
where, which are imposed upon the stiffener. The 
curve AB shown in the figure represents the mean de- 
flection of the wrinkled skin, assumed to be the same as 
the deflection of the stiffener. Assuming that the length 
of the stiffener is sufficient so that the skin will act as a 
membrane over the distance AB, it follows that there 
must be certain distributed loads lV’,., as shown, such 
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Fic. 23. Stiffener axial load failure. 
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Fic. 25 (left). Curved panel frame loads. Fic. 26 (right). Stif- 


fener eccentric bending moment. 


that the combined bending moment of W, and P,, about 
any point of AB is zero. Furthermore, the resultant at 
any point must be tangent to AB and equal to P;. In 
the"free-body diagram of the stiffener, CD is the cen- 
troidal axis. In this diagram also, the resultant load at 
any section of the stiffener is tangent to AB and equal to 
P,. Consequently the moment about any point of CD 
due to the load P, is equal to, but no greater than, Prz. 


The column effect is therefore negligible. 
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Stiffener local bending loads 


In the case where the stiffeners are far enough apart to 
invalidate partially the above argument, because the 
web no longer bends with the stiffener near the panel 
center, the transverse tensile stress fr begins to acquire a 
component tending to oppose lateral deflection so that 
the possibility of column failure is again inhibited. 

The effect of the tension field in partially supporting 
the stiffener was first pointed out by Wagner,’ who ac- 
counted for it by using a reduced column length in com- 
puting the strength of the upright. 

Open section stiffeners on beam webs show a tendency 
to fail by torsional instability. Whether or not this 
type of failure is opposed by action of the web is not 
clear. 

(4) Stiffener Secondary Bending.—The longitudinal 
and transverse tensile stresses f, and fr produce 
secondary bending stresses in the flanges and end 
stiffeners of beams. These stresses are of particular im- 
portance in the flanges of beams having widely spaced 
stiffeners. 

(5) Local Bending Loads.—This is the term here 
used to denote the bending loads imposed on the stringer 
by the buckling of the web. The effect is illustrated in 
Fig. 28. The loads per inch Q, produced on the stiffener 
are localized and sinusoidal in character. The maxima 
occur at points where the wrinkle crests intersect the 





ri’ 
ad 
at 
fr 


m 


fai 
ri\ 


St 


mi 
co 





~ STIFFENER 


yart to 
se the 
panel 
juire a 
) that 


orting 
10 ac- 
com- 


lency 
t this 
s not 


dinal 
duce 

end 
r im- 
aced 


here 

inger 
ad in 
ener 
Kima 
the 


fr. In the panel shown in Fig. 





DESIGN OF STIFF 
bending moments produced in the 
The loading shown 
in Fig. 28 is due to one panel only. If there are two ad- 
jacent panels, then the total stiffener load is the sum of 
the loads contributed by both panels. The resulting 
load on the stiffener is designated W,. The stiffener 
moments are larger if the two loads are in phase than if 
In particular, if the two ad- 


The 


stiffener. 
stiffener are in phase with the load. 


they are out of phase. 
iacent panels are identical and equally stressed, then the 
stiffener moments are maximum if the loads are in 
phase and identically zero if they are’ 180° out of phase. 
Fig. 29 shows a stiffener failure apparently due largely 
to local bending moments induced by the buckled web. 

(6) General Bending Loads.—These are the loads im- 
posed on the longerons or curved fuselage panels as a re- 
sult of the curvature of the tension field and are, in 
fact, identical with the loads W,; shown in Fig. 25. Re- 
actions are provided by the frames as shown in Fig. 30. 
The stringer is thus a uniformly loaded continuous 
beam, with the maximum bending moments occurring at 
Fig. 31 shows a curved panel in which the 
the 


the frames. 
longitudinal stiffeners failed under the action of 
general bending load. 

) Torsional Loads.—The local loads imposed on a 
stringer by two adjacent panels may combine to produce 
torsion instead of bending. This effect is particularly 
pronounced where the panels are identical and equally 
stressed, provided they buckle 180° out of phase as in 
Fig. 32. Fig. 32 also shows the resulting torsional loads 
on the stiffener. 


Rivet Failure 

The following loads contribute to rivet failure: 

(1) Shear Loads.—The total shear on a transverse 
riveted seam or attachment is found by vectorially 
adding the stresses f, and f,. On a longitudinal seam or 
attachment, the total stress is the vector sum of f, and 
33, the riveted attach- 
ment along the lower (transverse) edge failed in shear. 

(2) Tensile Loads.—The same loads that produce 
stringer bending also tend to cause failure of skin-to- 
stringer rivets in tension (see Fig. 34). Actually, the 
general load IV; always acts radially inward; conse- 
quently, it always reduces the rivet tensile load. It is 
the local load W, which is primarily responsible for 
panel failure by rivet tension. The total rivet load is 
equal to (WW, — W,)Spr. However, it should be noted 
that, in curved panel tests, the panel frequently flattens 
out, thus reducing the load Wg, and causing rivet tensile 
failure to occur under the action of W,. Fig. 35 shows a 
rivet tension failure. 


Stress Computation 

It is believed that the various effects discussed above 
can be adequately evaluated on the basis of the present 
method. No example is given in what follows of the 
computation of loads in curved panels, but the flat panel 
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Fic. 29. Stiffener local bending failure 
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Fic. 30. Curved panel stiffener general bending loads 





Curved panel stiffener general bending failure 
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Stiffener torsion loads. 





Rivet shear failure. 


Fic. 33. 


example should serve to indicate the type of calculation 
required in curved panel analysis. 


Analysis of Flat Panels 


In flat panel analysis, the basic geometric data and 
loading are determined, values of the panel stresses are 
read from the charts for R, = 1 and Rr = 1, these 
readings are multiplied by the correction factors corre- 
sponding to the true stiffening ratios, and the final 
stresses and loads are obtained by simple calculation. 


Panel Dimensions 
The panel 
Similarly, 


Fig. 36 shows the basic geometric data. 
length L is the length of sheet free to buckle. 
W is the buckled width. The longer dimension is al- 
ways denoted L. In the case of single web stiffeners 
with a single row of rivets, the web buckles all the way 
to the rivet line, and W (or L) is measured from rivet 
row to rivet row. If the web stiffeners are double, the 
strip of skin between stiffeners should be considered un- 
buckled, and IV is measured to the stiffener edges. 


Shear Stress 
In the case of a beam, the shear stress is taken as the 
total load divided by the distance between the cen- 


APRIL, 1950 
troids of the upper and lower riveted attachments time 


the sheet thickness. 


Effective Modulus 

Fig. 37 gives the effective web modulus of elasticity a. 
a function of the shear stress for 24STAL and 75ST 4] 
materials. 
flection data derived from a large number of shear panel 


This empirical relation is based on load-de 
tests made at Douglas Aircraft. The use of thes 
curves to account for panel behavior above the elast 
limit of the web is justified by the agreement found iy 
the previous section between test and analysis for th 
two kinds of material. 

It is logical also to use some sort of effective modulys 
for the stiffening elements, but this would require a solu. 
stiffener 
However, such a 


tion by successive approximations, since 
stresses are not known in advance. 
proc-dure can normally be avoided, since stiffeners, 
being in compression, are usually not so highly stressed 


as webs. 
Stiffening Ratios 
Basically, the stiffening ratios are given by 


- A, EK, a Ar Kr 


R = 
-” Wt E* , Li E 











it 4 Lge hy = We 






































(W,-W,) Sp 


Rivet tension loads 
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In these expressions, ‘A, is the amount of longitudinal 
stiffening associated with the area IVt of a transverse 
web section. Similarly, A; is the transverse stiffener 
area associated with the area Lt of a longitudinal web 
section. 

In the case of an interior panel of a beam (e.g., 
not the end panel), A, is the area of one stiffener where 
there are single stiffeners or is the area of two stiffeners 
where there is double stiffening (see section BB, Fig. 
36). Aris the total area of both flanges. 

In determining the longitudinal stiffening ratio for a 
beam, the bending flexibility of the flanges should be 
accounted for. This may be done by modifying the 
expression for R, as follows: 


(A, Wt) (Ez E 
WA, ( l ] ) Ey 
720L Ir, Ir, ky 


where /7,, and J;, are the moments of inertia of the 


R, 


upper and lower flanges. This equation was derived by 
finding the area of an effective stiffener which, when 
combined with infinitely rigid flanges, gives the same 
average strain as the actual stiffener-flange combination 
under the same loading. 

Where there are stiffeners on one side of the sheet 
only, the longitudinal stiffener area should be replaced 


by an effective value given by 


A, 


A;, = 
[1+ (Z/p)?) 


where Z is the distance from the stiffener centroid to the 
centerline of the web (see Fig. 36). This expression is 
due to Kuhn.? A similar correction may also be applied 
to transverse stiffeners. Where the stiffeners are 
double, no correction of this kind is required, but A, 
Should include the area of the strip of web between 
stiffeners. 

Note that the value of / used in the stiffening ratios 
should be the effective value of the web modulus. 
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Example of Flat Panel Analysis 


Determine the stress distribution in the beam shown 
schematically in Fig. 38. The design moment ts 14,500 
in.lbs. (compression on top), and the design shear is 


1,225 Ibs. per in. 
f, = 1,225/0.051 = 24,000 Ibs. per sq.in. 


From Fig. 37, 


E/E 0.487 
E = OAST XK 10.6 K 10° = 5.16 XK 10° Ibs. per 
sq.in. 
f. Wr 24,000 KX 8 XK 13.5 —_ 
Et? 5.16 X 10°X 0.0512 
L/W = 13.5/8 = 1.69 
0.127 rn 
A, = = (0.0558 sq.in. 
1 + (0.423 /0.363)? 
1 fh 0.538 KX 10.6 X 10° 0.971 
WE 8X0051*5.16X108 








0 
fe) 5 iO 5 20 25 30 35 
(f, 1000 PS) 
Fic. 37 Effective web modulus of elasticity 
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Fic. 38. Beam example 
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8° X 0.0538 
720 X 13.5 


al 1 
: -i + 1 = 1.130 
= = ae 
R, = 0.271/1.130 = 0.240 
0.281 + 0.375) X 10.6 X 10° ; 
Rr = ia area ——— = —— = 1.960 
13.5 X 0.051. X< 5.16 X 16 <X 10° 
(See Table 1.) 
. TABLE 1 : 
4 Uncorrected Corrected 
Type of Stress Value Value 
(o/fe) — 1 0.417 0.910 0.380 
L/fs 0.415 0.515 0.214 
T ‘Ss . y 0.422 1.235 0.522 
ae G)” (FF) : 26.5 1.015 X 1.025 27.6 
7 (*) (ey) ‘ 3.63 0.985 X 0.905 3.24 
SWLt \P./ 2 


Web 

a, = (0.880 + 1)24,000 = 33,100 Ibs. per sq.in. 
This figure should be multiplied by a factor to account 
for stress concentrations at the riveted attachment. A 
suggested correction is 


mg taa) = 
= GO _ > 
‘\s.-V2D 


Fs sign 
where Sz = rivet spacing, D = rivet diameter, and 
K = concentration factor = 1.1 for round-head 
rivets. 
Stiffeners 
P, = fitwW = 2,100 lbs. 
MM, = 3.24 EWLt/(E/f.)'“(WL/t?)"* = 40.6 in.Ib. 
E/fs = 5.16 X 10°/24,000 = 215 
WL/t? = 8 X 13.5/(0.051)? = 41,500 
A, = 127 sq.in., I, = 0.0167 in.4 
Z = 0.398 + (0.051/2) = 0.423 in. 
C = (0.398 or 0.602 


Stress in leg attached to web: 


2,100 40.60.3898 2,100 X 0.423 X 0.398 
-— = 
0.127 0.0167 0.0167 
= —38,710 lbs. per sq.in. (compression ) 


fy = 


Caps 
Compressive load due to transverse stress: 


Pr = ftL = 8,630 lbs. 
P,/cap = 8,630/2 = 4,315 Ibs. 


Load due to design moment: 
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Py = 13,500/14.351 = 940 Ibs. (compression top cap 


The upper cap is critical. 
Compressive stress: 
fe = (4.315 + 940) /0.281 = 18,700 Ibs. per sq.in. 


Bending stress at stiffener: 


fo = F, WC/12I 

F, = 24,000 X 0.214 X 0.051 = 262 Ibs. per in, 

: 262 (8)? 0.903 sa eta 7 

Jo ssi 12 X 0.041 = dU; IS. per sq.in. (com. 
pression) 

fax = 18,700 + 30,800 = 49,500 Ibs. per sq.in 


(compression) 


Web-Stiffener Attachment 
W, = 27.6Et/(E/f,)3/4(WL/t?)3/, = 44.2 Ibs. per in 
tension /rivet = 44.2 X 0.75 = 33.2 Ibs. per rivet 


Web-Cap Attachment 
/ 0.214 X 24,000 
Srr = fe + fr? = 24,600 Ibs. per sq.in. 
Frr = 24,600 X 0.051 
‘ rivet load = (1,252/2) 0.75 = 470 Ibs. per rivet 


I| 


= 5,140 Ibs. per sq.in. 


= 1,252 Ibs. per in. 


Stiffener-Cap Attachment 


P, = 2,100 We. 


If the load P, is assumed to be transmitted entirely 
through the two bolts, then the load per bolt is 1,05 
Ibs. However, it appears to be good practice to con- 
sider that one of the stiffener rivets just below the ca 
helps to transmit the load, on the assumption that a 
small piece of web in the corner of the panel acts asa 
gusset. 

Note in the above example that the depth of panel 
L, was taken as the distance between the centroids oi 
the upper and lower riveted attachments, assuming 
that these marked the limits of web buckling. In Tabk 
| the correction factor for Wy, is given by Ky, = 


Kwr,K WLo- 


CONCLUSION 


A graphical method of applying the previously de: 
veloped shear panel equations has been presented, anda 
comparison has been made between experimentally de- 
Reason- 


able agreement was found, and, where discrepancies 


termined stiffener stresses and chart values. 


appeared, it was possible that extraneous influences 
may have affected the test data. Since all of the aspects 
of panel behavior treated by the present method att 
part of a unified mathematical structure, it seems 
reasonable to suppose that verification of the predicted 
stiffener stress is at least a partial verification of the r 
mainder of the analysis. 
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DESIGN OF STIFFI 


There is need for additional experimental investiga- 
tion with regard to both the distribution of stress and 
the determination of allowable loads for the various 
panel components. In particular, the bending loads 
pre duced on the stiffener by the buckled sheet should be 
measured and compared with analysis values. In this 
connection, it should be mentioned that the charts for 
W, and MW, give loads acting at the center of the 
stiffener, assuming conservatively that the crest of a web 
buckle intersects the stiffener at its midpoint. The 
same is true for the quantities W, and AV/,. 

Besides the loads W, and Wz, there are torsional 
stresses in the web itself, and these may have some 
effect on the stiffener. No explicit expression was given 
in reference 1 for web torsion, although one could be 
derived, and the energy of warping was accounted for 
in the analysis. 

Although the loads W, and WZ, were small in the ex- 
ample analysis, their importance should not be dis- 
counted, because they increase rapidly with sheet thick- 
ness, for constant panel dimensions and shear stress. 

It is concluded that a practical method of shear panel 
analysis has been developed. This method has a 
rational mathematical basis, is readily and _ easily 
applicable to engineering problems, and is capable of 
giving satisfactorily accurate results through the plastic 
range to failure of the panel. Although only the stress 
analysis of panels has actually been demonstrated, it is 
anticipated that additional investigation and _ refine- 
ment will establish the validity of the method as a 
means of predicting ultimate loads. It is further ex- 
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pected that the method will prove to be applicable not 
only to thin web beams but also to curved wing and 
fuselage panels. Some degree of success has, in fact, al- 
ready been achieved in predicting ultimate loads for 
curved fuselage panels, although the results were not 
in a form suitable for publication at the time of prepara- 
tion of the paper. 
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Midwestern Conference on Fluid Dynamics 


A Midwestern Conference on Fluid Dynamics, together with a meeting of the American Physics Society, Fluid 
Dynamics Division, is announced by the Department of Mechanical Engineering of the University of Illinois. This 
conference is to be held on Friday and Saturday, May 12 and 13, 1950, concurrently with the dedication of a new Me- 


The Institute of the Aeronautical Sciences is one of the sponsoring societies, and its members are cordially invited 


Sessions on the following subjects are being planned: Aerodynamics—External Flow; Mechanical Engineering— 
Internal Flow; Symposium on Hypersonic and Low-Density Flow Techniques; Turbulence and Aero-Physics; 
Astro-Physics, Computing Machines, Multiple Theory of Fluid Dynamics, and Meteorology; Chemical Engineering; 
and Hydrodynamics and Continuation of External and Internal Flow. 

Reservations for accommodations and for the various banquets must be made before May 1. Reservation forms 
may be obtained from Prof. J. R. Carroll, Jr., Mechanical Engineering Building, University of Illinois, Urbana, IIl 

















Dynamic Stability at High Speeds from 
Unsteady Flow Theory’ 


I. C. STATLERT 
Cornell Aeronautical Laboratory, Inc. 


ABSTRACT 


A study is made of the longitudinal dynamic response char- 
acteristics of a typical fighter-type jet-propelled airplane based 
upon unsteady flow theory. The 14 longitudinal dynamic sta- 
bility derivatives are evaluated as complex functions of reduced 
frequency and Mach Number and used to determine the fre- 
quency response curves of normal acceleration and pitching ve- 
locity. The results are compared to the corresponding responses 
computed on the basis of the ‘“‘quasi-steady”’ theory, which neg- 
lects the effects of frequency and the higher order derivatives and 
includes the classical ‘‘time-lag’’ of downwash in the form of an 
M; term. It is shown that the discrepancies between the two 
can be accounted for by using the actual phase lag of the down- 
wash, as determined by nonstationary flow theory, instead of the 
“time-lag.’’ However, the differences between the simple and 
the exact results are so small that the inaccuracies of using 7; are 
probably no greater than those of the original theoretical as- 
sumptions or the experimental data. 


INTRODUCTION 


7 PROBLEM OF DYNAMIC STABILITY has been greatly 
simplified, in the past, by assuming that the sig- 
nificant stability derivatives were independent of the 
frequency of oscillations. This assumption, plus the 
condition that the speed remains constant, reduced the 
solution of the longitudinal equations of motion of the 
airplane (except for phugoid motion) to a simple linear 
second-order differential equation with constant co- 
efficients. However, a flight-test method for deter- 
mining the dynamic stability of airplanes, which is 
being developed by the Cornell Aeronautical Labora- 
tory, Flight Research Department, under a contract 
with the Aerodynamics Branch, Aircraft Laboratory, 
Wright-Patterson Air Force Base, involves sinusoidal 
oscillations of the airplane and its control surfaces or 
an equivalent analysis of the transient response. AI- 
though the frequencies of these oscillations are much 
lower than flutter frequencies, the question remains 
whether they are high enough to require that the non- 
uniform motion of the airplane be considered. Hence, 
in an endeavor to answer this question, an analysis of 
the longitudinal response characteristics of an airplane 
to a sinusoidal oscillation of its elevators is made which 
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depends upon a knowledge of the aerodynamic forces 
acting on an oscillating surface. 

The analysis in this paper is based upon a modifica 
tion of Dietze’s' ? solution for the lift distribution over 
a harmonically oscillating thin airfoil in subsonic com 
pressible flow. In addition, the effect of the oscillating 
wing wake upon the horizontal tail surfaces is con- 
sidered. The stability derivatives are found to by 
complex functions of frequency and Mach Number 
These results are then used to obtain the response 
curves for a recent U.S.A.F. fighter-type jet-propelled 
airplane and are compared to the results of the simpli- 
fied theory. A study of these results enables us to 
answer the following questions: 

(1) Are the effects of frequency in the range in which 
we are interested—i.e.,0 <w< 8 
sufficient magnitude to justify the additional compli- 


rad. per sec.—of 


cations of their inclusion ? 

(2) Are the higher order derivatives important? 
These include such terms as angular acceleration, trans 
lational acceleration, elevator velocity, and _ elevator 
acceleration, which are usually neglected in the simpli- 
fied theory. 

(3) Is the assumption of a “‘time-lag’’ term sulli- 
ciently accurate for determining the effect of the wing 
wake upon the horizontal tail ? 

The present paper is an abbreviated version of the 
final report published recently by the Cornell <Aero- 
nautical Laboratory under a similar title. This re- 
port presents not only the detailed derivations but also 
graphical representations of the results so that applica- 
tions can be made without an extensive knowledge of the 


theory. 


SYMBOLS AND CONVENTIONS 


Conventions 


Axes.—Stability axes are used throughout. These axes at 
fixed in the airplane and therefore move with it during a maneu 
ver. The origin 0 lies at the center of gravity of the airplane 
OX and OZ lie in the plane of symmetry, and OY is perpendicular 
toit. OX lies parallel to the free-stream wind vector (or parallel 


to the flight path) in trimmed level flight 


XY -axis—positive forward. 
)-axis—positive along the right wing 
Z-axis —positive downward. 


Linear Displacements.—A linear displacement along a positive 
reference axis is considered to be positive. 

Angular Displacements and Moments. 
ment or moment that is clockwise when viewed from the origi! 


An angular displace 


looking along a positive reference axis is considered to be pos! 


tive. 
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DYNAMIC STABILIT 


Accelerations.—Velocities and accelerations, 


and 
either linear or angular, are considered positive in the same sense 


Velocities 


.s the corresponding displacements. 
Control Surface Deflections.—Positive elevator angle is asso- 
ciated with a downward movement of the elevator trailing edge. 


Symbols 
U = true air speed, ft. per sec. 
w = incremental velocity along Z, ft. per sec. 
w = incremental acceleration along Z due to a change 
in w, ft. per sec.” 
q = pitching velocity, rad. per sec. 
q = pitching acceleration due to change in gq, rad. per 
sec.” 
a = elevator deflection, deg. or rad. 
j, = elevator velocity, di,/dt, deg. or rad. per sec. 
‘. = elevator acceleration, d*5,/dt?, deg. or rad. per 
sec.? 
D = differential operator = d/di 
g = acceleration of gravity, ft. per sec.” 
iz = moment of inertia about principal Y-axis, slug-ft.? 
m = mass of airplane, slugs 
S = wing area, sq.ft. 
c = mean aerodynamic chord, ft. 
l = length of fuselage, ft. 
w = angular frequency, rad. per sec. 
p = mass density of air, slugs per cu.ft. 
® = phase angle, general—phase relations indicated by 
subscripts—i.e., Pm, = phase of m with respect 
ton 
k = reduced frequency = wc/2U, nondimensional 
Fr = transverse force, lbs. 
r = fuselage radius, ft. 
N> = normal acceleration per g, nondimensional 
Z = force along Z, lbs. 
M = pitching moment, ft.lbs. 
1 OZ 1 OZ 
Zz = m Ze = » ee. 
m Ow m Og 
1 omM 10M 
M, - , == , ete. 
I, Ow Iy Oq 
dCi/da = —(2m/pUS)Z,, 
dC,/di, = —(2m/pU*S)Zs, 
dCn/da = (21,/pUSc)My 
dCm/db. = (2I,/pU?Sc) Ms, 


A bar placed over a variable indicates the maximum value of 


that variable during an oscillation. 


THE PROBLEM 


The problem, then, is to obtain theoretical expres- 
sions for the responses of an airplane in longitudinal 
sinusoidal motion at high subsonic speeds from un- 
steady flow theory. The total aerodynamic forces 
and moments are obtained by summing the effects of 
the five component parts of the overall problem. These 
are: (A) the oscillating wing, (B) the oscillating hori- 
zontal tail, (C) the horizontal tail in the wake of an 
oscillating wing, (D) the oscillating fuselage, and (E) the 
jet. Each of these problems will be discussed indi- 
vidually. 


A) The Oscillating Wing 

The first attempt to calculate the air forces on an 
oscillating airfoil in subsonic compressible flow was 
undertaken by Possio* and resulted in approximate 





bo 
w 
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Y AT HIGH SPEEDS 
values for the lift and moment on a wing without a 
control surface. Possio’s two-dimensional derivative 
theory applies for coupled simple-harmonic normal and 
pitching oscillations of an infinite, thin airfoil moving 
at subsonic speeds without shock waves. The flow is 
assumed to be adiabatic, and the Kutta-Joukowsky 
condition is satisfied at all times. From the linearized 
equation of motion for two-dimensional unsteady com- 
pressible flow and with the additional assumption of 
small disturbances, a simple differential equation is ob- 
tained in terms of an acceleration potential. The solu- 
tion is then represented by a special system of doublets 
with a chordwise distribution of intensity such that the 
correct normal velocity is obtained at all points on the 
profile. Since the lift distribution is directly pro- 
portional to the dipole strength distribution, the re- 
sultant integral equation also determines the lift dis- 
tribution. However, in order to determine this dis- 
tribution, it is necessary to resolve a singular integral 
equation with such a highly complex nucleus that the 
difficulties entailed are insurmountable. Possio’s ap- 
proximate solution is obtained by satisfying the bound- 
ary conditions only at a finite number of points along 
the airfoil. 

Another solution, devised by Dietze," ? is an iterative 
procedure based on the assumption that the lift dis- 
tribution is composed of the incompressible lift and a 
small additional lift caused by the influence of compres- 
sibility. Kiissner* has solved the integral equation 
for the incompressible case, and the additional term 
is obtained by repeatedly solving an integral equation 
with a kernal of the incompressible problem. In this 
manner, the lift distribution is obtained in terms of re- 
duced frequency, Mach Number, and the motion of the 
profile. The accuracy of this solution depends not only 
upon the number of steps for which the process is car- 
ried out but also upon the reduced frequency and Mach 
The convergence becomes worse with in- 
How- 


Number. 
creasing reduced frequency and Mach Number. 
ever, since we are interested only in extremely small 
values of reduced frequency, 0 < k < 0.1, and since the 
reduced frequency decreases with increasing Mach 
Number, for a given angular frequency, the lifts and 
moments for the oscillating wing with flap used in this 
paper are reasonably accurate. 

The two-dimensional results are corrected for three- 
dimensionality by a function corresponding to Reiss- 
ner’s® finite span correction factor in the incompressible 
case. This is justifiable since, in effect, we are correct- 
ing the sectional lift for compressibility and then modi- 
fying for finite aspect ratio as Reissner has done for the 
incompressible problem. 

The linearization of this problem makes it possible 
to divide complicated motions of the wing into so-called 
degrees of freedom from which the final solution is ob- 
tained by superposition. Accordingly, the simple case 
of harmonic oscillations of a rigid wing can be built up 


from: (a) translation of the wing as a whole, which 
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corresponds to flight along a straight path with chang- 
ing angle of attack; and (b) rotation of the wing about 
some reference axis, which corresponds to flight along a 
curved path at constant angle of attack. 


(B) The Oscillating Horizontal Tail 


The aerodynamic forces derived in the preceding 
section for an oscillating wing can also be applied to the 
horizontal tail after introducing the proper area, chord 
length, reduced frequencies, and efficiency factor. 
Since the longitudinal oscillation of the airplane is ob- 
tained through a sinusoidal motion of the elevators, 
it is also necessary to consider the forces and moments 
deflection, velocity, and acceleration 
In addition, it must be remem- 


due to the 
of the elevator. 
bered that the empennage moves in a nonuniform 
stream. 


(C) ‘The Horizontal Tail in the Wake of an Oscillating 

Wing 

During unsteady motions of the airplane, the wing 
leaves a sheet of vortices of varying strength in its wake. 
The velocity induced by these vortices has a pronounced 
effect on the direction of the flow of air in the vicinity 
of the tail. Influencing the downwash at any point in 
the wake are the bound vortices, the shed vortices, and 
the trailing vortices. Given the strength distribution 
for each of these, the induced velocity at any point can 
be determined. The incremental velocity due to a 
single vortex is found by the Biot-Savart rule, and the 
total downwash is obtained by integrating over the 
proper regions using the Cauchy residue theorem at the 
singular points. By this method, we find the downwash 
angle at the tail as a function of angle of attack; pitch- 
ing velocity; angular and translational accelerations; 
aspect ratio; and frequency of oscillation of the wing, 
of ‘tail length, and of Mach Number. The additional 
force on the tail due to the oscillating wing can then be 
obtained. 


(D) The Oscillating Fuselage 


Thus far we have considered the components of the 
airplane which contribute the major portion of the total 
aerodynamic forces and moments. In addition to the 
wing and empennage, however, we should determine the 
influence of the fuselage upon the dynamic stability. 

Consider that the fuselage, at any moment of the os- 
cillation, is in uniform circular flight. Munk*® found 
that, for an elongated airship turning with an angle of 
attack, a, the transverse force and moment in curvi- 
linear flight are the same as in straight flight with the 
same angle of attack a. The theoretical conclusion 
that the damping due to pitching is zero has been experi- 
mentally verified, for all practical purposes, in a series 
of tests conducted by Cotter.’? With this assumption, 
calculations of the stability derivatives for the body 
alone become simple. 


Laitone’ and Tsien® have found that, for an extremely 
slender body of revolution, the loading was unaffected 
by compressibility. ; 

Therefore, the transverse force distribution on a body 
of revolution moving at subsonic velocity is obtained 
by using the exact expression for the incompressible 
potential flow about a prolate spheroid: 


Fr f 4rar(dr/dx) 1 
(1/2)pU? ~~ Jo 1 + (dr/dx)? os 
(E) The Jet 


Since this method of obtaining dynamic stability char- 
acteristics from unsteady flow theory is applied to a 
typical jet-propelled airplane as an illustrative example, 
we should consider the effect of the jet. An aerojet, 
such as is used on this type of airplane, always damps 
the oscillations of an airplane, provided the distance 
of the air intake from the center of gravity of the air- 
plane is less than the distance of the exhaust from the 
center of gravity. Since the influence of the jet is de- 
rived from a Coriolis effect, the result is an increment in 
the moment due to pitching velocity. Using the re- 
lationships derived by Braun,'® it can be shown that 
this damping influence is negligible in computing the 
response characteristics in this example, although it may 
have considerable effect on faster airplanes. 

Thus, by summing the effects discussed in the pre- 
ceding sections, the 14 longitudinal dynamic stability 
derivatives of the obtained. These 
include the lifts and 
attack; pitching velocity; 
accelerations; and _ deflection, 
celeration of the elevator. The derivatives found by 
this method are complex functions of Mach Number 


airplane are 
moments due to angle of 
vertical and _ pitching 
velocity, and ac- 


and reduced frequency. 

With these stability derivatives determined and 
with the assumption that the forward velocity remains 
approximately constant, the basic dynamic longitudinal 
equations of motion may be written in operational form 
as is shown in Table 1. Solving simultaneously and 
substituting zw for the operator D, since simple har- 
monic oscillations are being considered, expressions 
for the pitching velocity and the normal acceleration 
per degree elevator deflection and for their phases are 
obtained. By introducing the values of the complex 
derivatives into these equations, we are able to plot the 
responses and their phase angles against the fre- 
quency of elevator oscillation for each condition of 
Mach Number, altitude, and center-of-gravity posi- 
tion. 

A comparison of the stability derivatives considered 
in the exact theory, the simple theory, and a modified 
simple theory, which is discussed below, is also pre- 
sented in Table 1. The exact theory contains all 14 
of the dynamic longitudinal stability derivatives as 
complex functions of frequency. 
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THE SIMPLIFIED THEORY 
In the past, the following assumptions were made 
which greatly simplified response calculations: 
(1) The stability derivatives are real and independ- 
ent of frequency. 
(2) The higher order aerodynamic 


Linke 23 Zi, 25, Mj, Mj,, and Mi, are negligible. 


derivatives, 


(3) The downwash lag can be approximated by the 
time required for the flow to travel from the wing to the 
tail. The downwash lag terms as derived by Cowley 
and Glauert!! depend on the translational acceleration 
dw/dt and are therefore contained in the acceleration 
derivatives \/; and Z,;. It is important to note, how- 
ever, that they arise from a different cause than that 
associated with the acceleration derivatives. 
although an J/, term is included in the simpli- 
fied theory, it is not comparable to the 1/; derivative 
obtained from unsteady flow theory, which is due solely 
to apparent mass effects resulting from the translational 


usually 
Hence, 


acceleration. 


With these simplifications, the equations of motion 
can be written in a form corresponding to a simple 
spring-mass system with damping, where the equivalent 
spring and damping constants are functions of the sta- 
bility derivatives. Laitone'? has shown that a knowl- 
edge of these constants allows the prediction of the 
airplane motion for any arbitrary control surface mo- 
tion. Furthermore, a simple evaluation of relative 
dynainic stability is provided by a comparison of the 
spring and damping constants for different airplanes. 


THE SUPERSONIC CASE 


The basic theory for stability and control at super- 
sonic speed is exactly the same as for the subsonic case 
i.e., the equations of motion as shown in Table 1 are 
valid if the correct values for the stability derivatives 
are introduced. Hence, the problem is to determine 
the 14 stability derivatives for an airplane that is pitch- 
ing sinusoidally at supersonic velocity. 

The evaluation of the forces and moments on an os- 
cillating airfoil in a supersonic stream has been ade- 
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quately covered by several authors'*~'® for the two- 
dimensional case. Since we are interested in extremely 
low values of frequency, a first-order approximation is 
valid and provides simple expressions. Furthermore, 
we are able to assume that the wing attitude changes 
slowly relative to the rate at which a disturbance is 
transferred downstream in the Mach cone. This per- 
mits the use of Evvard’s' method as an approximate 
three-dimensionality correction for the unsteady case. 

With the forces and moments on the wing known, 
those on the horizontal tail can be found in the same 
manner as in the subsonic case, except for the deter- 
mination of the effect of the oscillating wing wake. 
Theoretically, for a conventional airplane configuration, 
only that portion of the horizontal tail which lies within 
the wing Mach cone is affected by a downwash. This 
downwash has been determined by Lagerstrom and 
Graham" for the steady case. The same approach is 
utilized for the case of nonuniform motion. 

Quantitative results were not available at the time of 
writing; however, it can be stated that, to a first-order 
approximation, the stability derivatives are complex 
functions with real parts independent of and imaginary 
parts proportional to the frequency. 

This problem of dynamic stability at supersonic 
velocity is treated in more detail in the forthcoming 
Cornell Aeronautical Laboratory report. 


DISCUSSION OF RESULTS 


The method of determining dynamic stability at high 
speeds from unsteady flow theory is applied in an illus- 
trative example. The stability derivatives obtained in 
this manner are plotted against reduced frequency for 
several Mach Numbers in Figs. 1-14. The curves 
show the effects of frequency and Mach Number upon 
the longitudinal dynamic stability derivatives of a typi- 
cal jet-propelled fighter-type airplane. Wind-tunnel 
values for the static aerodynamic derivatives are cor- 
rected by the theoretical frequency effects. 

It would seem from these curves that the influence of 
frequency becomes increasingly more important with 
This problem, however, is 
It must 


increasing Mach Number. 
not so critical as it seems to be at first glance. 
be remembered that the reduced frequency is inversely 
proportional to the velocity. Therefore, for a given 
frequency, the reduced frequency diminishes by the 
ratio of velocities for increasing Mach Number. 

In Figs. 15-22, the response curves* of angular ve- 
locity and normal acceleration per degree elevator de- 
flection and their phases are plotted against frequency 
for Mach Numbers of 0.30, 0.50, 0.70, and 0.75. Com- 
parison of the results of simple theory with exact theory 

* The response curves have been faired into the zero frequency 
values, although they are not correct in the phugoid frequency 
range, which is well below 0.5 rad. per sec. At these low fre- 
quencies, the speed change is greatest, and its effect becomes 


extremely important. 
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shows excellent agreement in both phase and amplitude 
except at maximum amplitude or resonance frequency 
Since the two theories check at the high frequencies, thy 
contributions of nonstationary flow and the higher order 
derivatives, which are expected to be greatest at th, 
highest frequencies, cannot be important in the rang 
0< w< Srad. per sec. 
appear at resonance frequencies can only be explained 
by such a term as the downwash lag which has its major 
effect when the change and rate of change of lift ar 
greatest. This supposition was verified by further 
calculations based on a procedure referred to in the pre 


Therefore, the differences that 


ceding sections as the modified simple theory. Thy 
responses were recomputed using only those derivatives 
considered in the simple theory with the exception oj 
the downwash lag term. Instead of \/;, the complex 
form of de/da as obtained from unsteady flow theor 
was used. The agreement with the so-called exact 
results is excellent. This is shown by the comparison 
between the circled points and the solid lines in Figs, 
15 99 


CONCLUSIONS 


In the introduction of this paper, three questions 
in the problem of dynamic stability at high speeds 
were raised. Consider the first two: (1) Should the 
effects of frequency upon the stability derivatives be 
included and (2) should the higher order derivatives be 
included ? 

It has been shown that, at the highest frequencies 
investigated when these two effects would be maxi- 
mum, there is excellent agreement between the simple 
and exact theories. Therefore, the contributions due 
to both are not important for accurate theoretical re- 
sponses at subsonic speeds, either because these terms 
themselves are negligible or because their effects are 
compensating. 

An explanation of the differences between the simple 
and exact theories must, therefore, lie in the answer to 
the third question: 

(3) Is the simple ‘‘time-lag’’ assumption sufficiently 
accurate? In answering this question, it is important 
to consider the degree of accuracy desired. Undoubt- 
edly, the complex form of the downwash as obtained 
from oscillating airfoil theory is the more correct, since 
the transient effects of the wing lift are included. The 
responses obtained by the exact theory have been dupli- 
cated merely by using the complex downwash instead 
of the classical ‘‘time-lag’’ term, /;. However, the 
differences between the simple and the exact results are 
so small, even at the resonance frequencies, that the 
inaccuracies of using 1/7; are probably no greater than 
those of the original theoretical assumptions or the ex- 
perimental data. 

In conclusion, it can be stated that the quasi-steady 
method of computing response characteristics at sub- 
sonic speeds is reasonably accurate, at least for longi- 
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tudinal oscillations. If the more correct values of the 
downwash as obtained from unsteady theory are to be 
used instead of the “‘time-lag’’ /;, term, they can be 
obtained easily from general charts as functions of re- 
duced frequency, Mach Number, and tail position with 
respect to the wing, which will be contained in a forth- 
coming Cornell Aeronautical Laboratory _ report. 
Hence, the determination of the longitudinal dynamic 
stability of a conventional airplane in subsonic com- 
pressible flow is no more difficult than the equivalent 
incompressible problem. 

The possibility of unsteady flow effects upon the 
lateral dynamic stability characteristics of an airplane 
is still open to question. In particular, the problem 
of lateral “‘snaking’’ oscillations may be clarified by 
further investigation based on nonstationary flow 
theory. A preliminary study, however, reveals no 
obvious reasons for the author to suspect greater effects 
in the lateral case than were presented in the longi- 
tudinal case. 


REFERENCES 


! Dietze, The Air Forces on the Harmonically Vibrating Wing 
Compressible Medium at Subsonic Velocity (Plane Problem), 
AAF Translation F-TS-506-RE, November, 1946. 

2 Dietze, Part IT: Numerical Tables and Curves, AAF Trans- 
lation F-TS-948-RE, March, 1947. 

3 Possio, C., The Aerodynamic Forces on an Oscillating Profile 
in Compressible Fluid at Subsonic Speed, Aerotecnica, Vol. 18, 
No. 4, pp. 441-458, April, 1948 (R. T. P. Translation No. 987). 

‘ Kiissner, H. C., Zusammenfassender Bericht uber den in- 
stationaren Auftrieb von Flugeln (summary report on the non- 
stationary lift on wings), Luftfahrtforschung, Vol. 13, No. 12, 
pp. 410-424, December, 1936. 


APRIL, 1950 





sp 1 —_ 


CG. at 27% MAC. | — 
Altitude = 20,000 Ft 


y M = 
- J Di iach Number 075 | 


8 
NX 





















































fe) 
) 2 oS... 10 rom 
| —__ Exact ] 
Ya F Rese Simple 
S © Simple with a 
complex 9%, | | 








DEG. LAG 
é 
Z 








€ \ 


es 
































fe) 2 4 
a 8 10 


Normal acceleration per degree elevator deflection ys 
frequency. 


Fic. 22. 


® Reissner, E., Effect of Finite Span on the Airload Distribu- 
tions for Oscillating Wings. [—Aerodynamic Theory of Oscillating 
Wings of Finite Span, N.A.C.A. T.N. No. 1194, March, 1947. 

6 Munk, M., The Aerodynamic Forces on Airship Hulls, 
N.A.C.A. Report No. 184, 1924. 

7 Cotter, W. E., Jr., Summary and Analysis of Data on Damp- 
ing in Yaw and Pitch for a Number of Airplane Models, N.A.C.A. 
T.N. No. 1080, May, 1946. 

5 Laitone, E. V., The Linearized Subsonic and Supersonic Flow 
About Inclined Slender Bodies of Revolution, Journal of the Aero- 
nautical Sciences, Vol. 14, No. 11, pp. 631-642, November, 
1947. 

® Tsien, H., Supersonic Flow over an Inclined Body of Revolu- 
tion, Journal of the Aeronautical Sciences, Vol. 5, No. 12, pp. 
480-483, October, 1938. 

1 Braun, G. W., Damping Influence of an Aero-Jet Dud, 
ATSC-(T-2), IRE-22, December, 1945. 

1! Cowley, W. L., and Glauert, H., The Effect of Lag of the 
Downwash on the Longitudinal Stability of an Airplane, R. & M. 
No. 718, February, 1921. 

12 Laitone, E. V., A New Approach to the Analysis of Airplane 
Motions, Comptes Rendus du Sixieme Congress International de 
Mecanique Appliquée, September, 1946. 

13 Possio, C., The Aerodynamic Action on an Oscillating Aero- 


foil at Supersonic Speed, Acta Pontif. Acad. Sci., Vol. 1, No. 11, 


pp. 93-105, 1937 (A.R.C. Translation No. 7668). 

‘4 yon Borbely, S., Aerodynamic Forces on an Harmonically 
Oscillating Wing at Supersonic Velocity (2-Dimensional Case), 
Z.A.M.M., Vol. 22, No. 4, pp. 190-205, August, 1942 (R.T.P. 
Translation No. 2019). 

16 Garrick, I. E., and Rubinow, S. I., Theoretical Study of At- 
Forces on an Oscillating or Steady Thin Wing in a Supersoni 
Main Stream, N.A.C.A. T.N. No. 1383, July, 1947. 

6 Evvard, J. C., Distribution of Wave Drag and Lift in the 
Vicinity of Wing Tips at Supersonic Speeds, N.A.C.A. T.N. No 
1382, July, 1947. 


(Continued on page 2535) 





rr, 


Ex 
comy{ 
disco 
at a | 
tion 
show 
(n) @ 
ticul: 
funct 
tion ¢ 
solut 
alone 
Two: 
at th 
patte 
veloc 
and d 
soluti 
soluti 
desire 


sign 


T 


andr 
mom 
is do: 
and ¢ 
the 1 
neces 
An 
cade 
Cann 
ods b 
ties 1 
may 
dime: 
assoc 
veloc 
angul 
fluid 
to inf 
ever, 


ment 


Rec 


Ml Vs. 


tribu- 
lating 


lulls, 


amp- 
C.A. 


Fl ow 
Lero- 


iber, 


volu- 
pp. 


uct, 


” the 
, M. 


lane 
1 de 


ero- 


ll, 


the 
No. 








Some Exact Solutions 


of the Flow Through 


Annular Cascade Actuator Discs’ 


STEPHEN L. BRAGG? ann WILLIAM R. HAWTHORNE? 
Massachusetts Institute of Technology 


SUMMARY 


Exact solutions are presented for two types of flow of an in 
compressible inviscid fluid through an annular passage with a 
discontinuity in tangential velocity at one section, as would occur 
at a cascade of airfoil blades of negligible chord. When stagna 
tion conditions are constant throughout the flow region, it is 
shown that the product of the tangential vorticity component 
(n) and the radius (7) is constant along a stream line. The par- 
ticular case when this product is a linear multiple of the stream 
function (y) is solved by dividing y into persistent and perturba- 
tion components and separating the variables. 
a function of radius 


A second exact 
solution is obtained when the product (nr) is 
alone and when the stagnation conditions are allowed to vary. 
Two solutions of either of the above kinds can be exact!y matched 
at the discontinuity section, thus determining the complete flow 
pattern for arbitrarily selected conditions at infinity. The axial 
velocity at this section is almost exactly the mean of upstream 
and downstream values at the same radius. The family of exact 
solutions may be of use in estimating the accuracy of approximate 
solutions and, as shown in a numerical example, may match some 
desired flow conditions sufficiently closely to be of interest in de 


sign 
INTRODUCTION 


_ ESSENTIAL FEATURE of axial compressors and 
turbines is the arrangement of alternate stationary 
and rotating cascades of airfoils; these alter the angular 
momentum of the working fluid as it passes so that work 
is done on or by it, thus raising or lowering its pressure 
and energy levels. In designing such turbomachinery, 
the required work distribution is known, and it is 
necessary to find the blade profiles that will produce it. 

An analysis of the flow through a single annular cas 
cade is a first step in the study of this problem. This 
cannot be done by the classical two-dimensional meth- 
ods because of the radial pressure gradients and veloci- 
ties induced by the fluid rotation, although attempts 
may be made to correlate data with results from two- 
dimensional cascades. The radial pressure gradient is 
associated with centrifugal effects due to the tangential 
velocities produced by the blades. Conservation of 
angular momentum demands that, in the absence of 
fluid shear, these persist as the fluid flows downstream 
to infinity after a cascade. The radial velocities, how 
ever, are purely local phenomena caused by the adjust- 
ment in mass flow distribution associated with the dif 
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ferent velocity profiles on either side of the cascade. 
We may then regard the flow as divided into two parts 
upstream and downstream. In each, the flow condi- 
tions are basically those that would obtain at infinity, 
but, in the region of the cascade, perturbation compon- 
ents are introduced to smooth out the transition from 
one region to another. 

Various authors have suggested analytical approaches 
to the determination of the flow through cascades in 
order to obtain the blade profiles for a given duty. 
Meyer! has studied the flow through free vortex blades. 
Marble* has presented a general solution, but this in- 
volves some mathematical approximation in order to 
linearize the equation of motion. 

This paper describes two physically possible types of 
incompressible nonviscous flow for which exact solu- 
tions of the equations of motion are found and applied 
directly to the limiting case of a cascade of negligible 
chord width. This is a special case of the actuator disc 
of propeller theory. The flow is assumed to be steady 
and uninfluenced by external forces, except the nor- 
mal forces at the parallel walls of the annulus. At one 
particular section, the actuator disc, there is a dis 
continuity in the tangential—though not in the other 
components. Rotational 
there is no variation in 


velocity symmetry is as- 


sumed—that is, flow condi 
tions at a given radius with angular position. 
is, in effect, the limiting case of a cascade containing an 
infinite number of blades of infinitesimal chord which 


The disc 


exert an impulse on the fluid in the tangential direction. 
The basic equations employed are Euler's Equation of 
Motion, the equation of continuity, and the thermody- 
namic entropy and energy equation. Heat transfer 
to the fluid is neglected, and it is assumed that the 
fluid flows isentropically from a region where it is all at 


the same entropy. 


Basic EQUATIONS OF ANNULAR FLOW 


In the absence of body forces, Euler’s equations of 
motion for rotationally symmetric steady flow of a non 
viscous fluid referred to cylindrical polar coordinates 


r, 6, and z (Fig. 1) may be written: 


Ov Oc, Ca lop , 
( + ¢, -- = — (radially (1) 
oO? Oz } por , 
OC, Oly , CCe 
c c = () (tangentially (2 
or Oz } 7 
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fr, Cy, and c, represent the velocity components in the 
three coordinate directions. 
The equation of continuity under similar conditions 


becomes: 
(4) 


The components of vorticity are defined by the rela- 


(0/0z) (prc:) + (0/Or) (prc,) = O 


tions: 
— = —Oc,/0z (radial component) (5) 
n = (0c,/0z) — (Oc,/Or) (tangential component) (6) 
¢ = (1/r) (0/dr) (re,) (axial component) (7) 


The enthalpy of the fluid at a point, 4, and the stag- 
nation enthalpy, //—i.e., the value of the enthalpy that 


would be reached if the fluid were brought to rest adi- 


abatically—are connected by the steady flow energy 
equation: 
= h + (1/2) (c,? + cp? + c,”) 

If we assume that all the fluid was initially at the 
same entropy and that all subsequent changes occur 
isentropically, the thermodynamic relation 6h = 6p/p 
applies, and so, for all variation through the flow, 

611 = (6p/p) + (1/2)8(c?) (8S) 

Combining Eqs. (1), (2), and (3) with Eq. (8) to 


eliminate pressure differentials and introducing the 
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vorticity components from Eqs. (5), (6), and (7), 


§¢g — ne, = OH/Or (9) 

fc, — &, = 0 (10 

nC, — {Cg = OH/Oz (1] 

Eq. (4) shows that it is possible to define a “Stream 
Function,”’ y, of ry and z, such that 

prc, = po(Ow/Oz) (12) 

prc, = —po(Ov/Odr) (13 


where po is a constant reference density. 
From the tangential equation of motion, Eq. (2 


[c,(0/Or) + c¢(0/0z)](rc,) = O 


This agrees with the physical concept that in the ab. 
sence of external forces angular momentum is pre- 
served, since it shows that rc, is constant along a stream 
#6 is a function of y. Also, H is a 
165). The relation 
between @ and the vorticity components is analogous 


line. Thus, rc, = 
function of W (reference 3, article 
to that between y and the corresponding velocities, for 

ré = —06/dz (14 
+ = O6/Or (15 


From Eqs. (6), (12), and (13) 


diT/dy = (1/r*)[0(d6/dy) (po/p)rn] (16 
0 l 0 ra) O 

= Yas + Ye 1) 244%) 
p r p* Oz or 


where y,, ., etc., denote partial differentiation of y with 
respect to r and g, etc. 
It is of interest to note the similarity between Eq. 


(16) and the relation derived by Crocco? for flow behind 
an axially symmetric shock, which is 
dH ds I 
= a = » rn (16a 
dy dy rp 


where s is the entropy and 7 is the absolute tempera- 
ture. 
Hence, Eq. (16) may be extended to a more general 
one allowing for entropy variation: 
Tds 1/.d6 
- (6 


— 4° ) (16b 
—s 


dy dy 
In incompressible flow, the relations (16) and (17 


simplify to 


dH/dy = r?)[0(d0/dy) + nr] (18 

where 
= rr + v2 — (1/r)¥, (19 
Features of various types of flow may be distin- 


guished by examining Eq. (18). 
with H = 


plies constant stagnation pressure 


Constant.—This im- 
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FLOW THROUGH 


flows through stationary or rotating cascades 


isentropic 


with constant circulation along the blades. 

Since dH /dy = 0, nr must be a function of y and con- 
stant along a stream surface. Also, from Eqs. (9), (10), 
and (11), 

(20) 


E/C, = /Cg = §/C, = Q/C 


where Q is the resultant vorticity. The vortex lines 
thus follow the stream lines. 

(b) Isentropic Flows with rc, Zero or Constant. 
From Eq. (15), 


dH/dy = n/r (21) 


and, hence, 7/7 is constant along the stream surface. 
The condition that »/r should be constant along a 
stream surface, mentioned by Marble? and in Lamb 
(reference 3, article 165), was applied to axially sym- 
metric flow— that is, flow with no tangential velocity. 
As shown here, it also results when rc, is constant and, 
hence, when £ and ¢ are zero. It ranges from the 
trivial, vortex-free case when 7 = 0 and J/7 = constant 
(the so-called free vortex flow) to flow with varying [7 


and 7. 


THE First Exact SOLUTION 
In order to obtain the first solution, incompressible 
fow with constant stagnation enthalpy has been as- 
Then, from Eq. (18), 
mr = fi) (22) 


The perturbation concept of the flow suggests that a 
suitable form for the stream function is given by the 


sumed. 


relation: 


yv=yviot+ De**™y,’ (23) 
n 
where Y, and y,’ are functions of 7 alone. &, is as- 
sumed positive, and the sign of the index +k, is such 
that the terms in the summation all vanish at great 
distances from the origin. 
Substituting in Eq. (19), 


—1/r¥. + De**™x 


(yp) = nr = Va 


l 
vn’) (24) 
r: 


», the summation terms vanish, leaving 


(e,"v a T aad — 


As |z 


IW) = Yor 


In the particular case where /(y) is a linear function 
of Y, it is possible to separate the variables by subtract- 
ing Eq. (25) from Eq. (24) and so solving for y,,’ and 
¥.. The case where f(y) is a constant will be considered 
later as part of the second exact solution. 

In order to nondimensionalize the equations, the 
teference length r,, equal to the annulus outer radius, 


and the reference velocity c,, are introduced. ,, is de- 
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fined as the mean axial velocity of the flow, which is 
constant along the annulus in incompressible flow, and 
is given by 


tip tip 
Ci 2arr dr = 2rrc, dr 
hub hub 


Cn = —2[W)inn (7,2 — r,”) 


whence, 


(26) 


Then, putting f(¥) = —w*(W/r,*), where uw is a con- 


stant, Eq. (25) becomes: 


Von > (l/r), + (u?/r?)W. = 0 (27) 
This has as solution a first-order Bessel function: 
Vo/PCm = —R(a/u)Ai(uR) (28) 


where R = r/r,and 

A,(uR) = J,(uR) + (a’ a) Y,(uR) 
a and a’ are constants, but only one is independently 
variable, for, eliminating c,, from Eqs. (26) and (28), 


a = p(l — R,?)/[Ai(u) — R,A\(uR,)} (29) 


where 7;,/7, = R,. 
For the perturbation component, subtracting Eq. 
(25) from Eq. (24) and substituting for f(y), 
+ knz , I , 9 we , 
de 7 Pics a Vnr + + 9 Vn an 0 
r r,? 
Omitting subscripts for clarity, it is seen that each 
term of the summation must satisfy the relation: 
Ver’ — (1/r)br’ + [R? + (u?/r)]y’ = 0 
Solving and putting 
2 _ k*r ,? aa Ti (30) 


another first-order Bessel function expression is ob- 


tained: 


W'/r?Cm = —R(b/dr)B,(AR) (31) 


where 
B,(AR) = J,(AR) + (6'/b) Y,(AR) (32) 


and } and b’ are constants of integration. The full 
solution for the stream function is, therefore, from Eq. 


(23), 


, ii 


r7Cm bl 


+h Dn 
— Doe***— RB(A,R) 333) 
n n 


Differentiating and applying Eqs. (12) and (13), 


C:/Cm = AAo(pR) + p e* ** bh Bo(X,R) (34) 


n 
- - - = i 
C)/Cmn = * é 


"(R11 An )O,Bol(d,R) (35) 


From Eqs. (18) and (27), 


6(d0/dy) = —nr = (p?/r/7)¥ 
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Integrating, 
#? = (u?/r 7) yp? oh Ti a 


where g is another nondimension?! constant of integra- 


tion. Therefore, 


a\? (eViv\, ge _ 
e i (‘) (+) i R? 
g 


b, +bpg ‘ 
| oAs(uk) +> % es Bidek) | + R? (36) 


MATCHING CONDITIONS AT INFINITY 


Putting the perturbation components equal to zero in 
the Eqs. (33) to (36) for the complete solution derived 
in the last section, the following conditions at infinity 
are derived: 


Cz2/€m sand aAy(pR) (37) 
[(u/R)/(Wo/ri?€m) |? + (g/R?) 
= —(a/p)RA,(uR) 


(Cog "se = 
Wiel fe Cm 


Three of the constants involved are arbitrary, the 


(3S) 


fourth being fixed by Eq. (29). It should, therefore, 
be possible to make a good approximation to any de- 
sired velocity distribution at infinity. 

A particular case is that of free vortex flow for which 
c, and rc, are constant—1.e., independent of radius. 


This condition is satisfied when a = 1 anda’ = yp = 0, 


and so X,, = k,r;. The full solution then becomes: 
C,/Cm =e 1+ > 0 DBA) (39) 
n 
C6 Joa = > e= 8/7) 5 Bi(X,R) (40) 
Nn 
G/ta=e/K* or Gc)* = 2(r2,.)* (41) 


: Toillustrate the possibility of matching Eqs. (37) and 
(38) to a prescribed distribution, a flow in which c, var 
ies as the square of the radius has been used, with c, = ¢ 
The 


exact equations for this flow at infinity are, denoting 


at the tip and a hub to tip radius ratio of 0.45. 


tip axial velocity by c,,: 


Gift = (42) 
Ce/t = V (3/2)[(5/3) — R‘] (45 

y r 2 /3 a ta l na oO 
rite | WING” Vag 4  N\3 R') (4 


where, from Eq. (26), 


Cn s. =e ( Fre a Va | 
Cae l — (0.45)*\7 Ce, 17Ce: 


These are plotted against radius ratio as the broken 


( 15) 


lines in Fig. 2. 
Values given by Eqs. (37) and (38) have been fitted 
to this curve by trial and error. The solid curve shows 


the result of taking values for the constants of 
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There is a slight curl in the rc, curve at the inner ra 
dius, but agreement is otherwise close, even in such an 


extreme case as the one chosen. 


MATCHING CONDITIONS AT OTHER BOUNDARIES 


The constants 6’,b and X are determined by the 
boundary conditions at the annular walls, which aré 
0 for R = 
This condition will be 


assumed parallel. 
R, and R = 


satisfied if, from Eq. (35), 


It is required that « 
| for all values of sz. 


B,(X,) B,(X,R, = () 


for all values of X,,.. Therefore, from Eq. (32), 
a” Ji(X») 
b, V(X, 


Ji (A, Rr) 
Vi(A,R, 
This provides an infinite number of possible values 


for 6’/b and X. 
Jahnke and Emde’s ‘‘Tables of Functions,’ 


The first six roots may be found in 
‘and highei 
roots, if they were required, could be worked out from 
the approximate formulas: 


Ji(x) = vWf 2/xx cos [x (32/4 


Vi(x) = VY 2/mx sin [x — (37/4 
whence, 


tan [A, — (3m/4)] = tan [A,R, — (39/4 
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FLOW THROUGH ANNULAR CASCADE ACTUATOR DISCS 


There fore, 


d, = Nw/(1 — R;) (47) 


where V is any integer. Note that, even for a hub: tip 
radius ratio as small as 0.5, the lowest value of 2, is 
about 6. The corresponding values of k, are obtained 
from Eq. (30). 

Only the values of the constants 6 are now unknown, 
and these depend on matching at the blade rows. If 
upstream and downstream flows are both defined by 
lutions of our form, exact matching may be obtained 
in the limiting case of an actuator disc. At the disc, 
continuity in radial and axial velocities is required at all 
radii, but there may be discontinuities in tangential 
components. If the disc is situated for convenience at 
: = 0, exact matching of the c, values will be obtained 
from Eq. (35) if 

Rnana = —Ruy Dry = Ron (48) 
where suffixes u and d denote upstream and downstream 
conditions, respectively. 

Matching of the axial velocities is obtained from 
Eqs. (34) and (37) if 


k, n k bn . - ~ 
> ( + ) By(X,,R) = (: ‘) —_ (‘ ‘) (49) 
n k nd k nu Cm u“ Cm d 


It is now possible to find &, and the other coefficients 
by a method analogous to that employed in determin- 
ing those of a Fourier series. Multiplying both sides of 
Eq. (49) by RB(A,,.R)dR and integrating from hub to 
tip, we find that all terms on the left-hand side for which 
m#nare zero. Performing the integration for m = n 
and making substitutions for c,,, from Eq. (37), we ob- 
tain the following expression for k,,,: 


] l l 
sho; +5 - yr [Bo?(A.,) = R,7By?(A,Rr) | aos 


nd 


(daba Ria?) [A i(ua) Bold») = R,Ai(uaR,) Bi(dA,R,) | - 
(Cyby Rea’? [A 1(u,) Bo(An) =— Ry (u,kR,) By (A,R,) | (50) 


This fixes all the coefficients in the series expansion, 
and thus the flow with given conditions at infinity up- 
and downstream is fully determined everywhere. 

One case is of special interest. If u, = wy, then 
Ree = bag, —Ony = Daa, anid Eq. (34) reduces to 


C; - 7 
| ‘) = (° >) + » by, aBo(\,R) - (‘ *) = 
Cm /a Cm 


m/ disc n u 


l * _ 
> bnaBo(A,R) = | (2) 4 (22) | (51) 
2 = Em d Cm u 


that is, the axial velocity at the disc is exactly the 
arithmetic mean of those at the same radius upstream 
and downstream at infinity. 

If u* is small compared with ),°, it will be extremely 
small compared with higher order \,, values. Thus, 


Te ae? 


bo 
_ 
“ 


AXIAL 
VELOCITY 


TANGENTIAL 
VELOCITY 


RADIAL 
VELOCITY 





1.0 0.9 08 O7 0.6 as 
RADIUS RATIO, + 
t 


Fic. 3. Conditions near the actuator disc. 


which again gives the result of Eq. (51). It seems 
likely that this relation will hold in many practical 
cases, so that, in general, the axial velocity at the disc 
will be almost exactly the arithmetic mean of the in- 
finity values at the same radius. 

On the other hand, y = y,, at root and tip radii and 
is little different elsewhere due to the 1/X,, factor in the 
coefficients of the series expansion [Eq. (33)]. Then 
Eq. (36) shows that there will be only a small differ- 
ence at a given radius between the tangential velocity 
just downstream of the actuator disc and the corre- 
sponding infinity value. 

To illustrate the application of the method, a flow that 
is axial and vortex-free upstream at infinity has been 
matched at the actuator disc to the approximate Cc, 
varying as r* example downstream at infinity already 
mentioned (shown as solid curves in Fig. 2). The first 
six roots of Eq. (46) for R, = rp/r, = 0.45 are: 


Nis 5.843 11.495 17.22 22.90 28.6 34.3 
bn’/bn —2.31 3.86 0.89 0.126 —0.488 —1.69 
and these give, with Eqs. (48) and (50), 
boy 0.0818 —0.0166 —0.0230 —0.0265 —0.0155 —0.0070 
bn,g 0.0856 —0.0167 —0.0230 —0.0265 —0.0155 —Q0.0070 


The calculated values of c,/c,, and c./c,, at the actu- 
ator disc and c,/c,, just downstream are shown in Fig. 3. 
The values of c./c, at the disc are almost exactly the 
arithmetic mean of upstream and downstream infinity 
values. Because of the large difference between the 
axial velocity distributions at infinity on either side 








248 JOURNAL OF THE AERONAUTICAL SCIENCES—APRIL, 1950 


of the disc, the radial velocities are high—up to 8 per 
cent of the mean axial velocity. Since rc, is constant 
along a stream line and the radial velocities are in- 
flows, the tangential velocities at the disc are less than 
those at the same radii at infinity downstream, except 
at the hub and tip where, the radial velocities being 
zero, no change occurs in Cy/Cm. 

The angle of flow, a, is also plotted. The effects 
mentioned above tend, in this case, to increase the angle 
at the tip and reduce it at the root as the flow moves 
downstream. An approximation suggested for rough 
design work was to calculate the flow angles near the 
disc using the mean axial velocity of upstream and 
downstream flows and the downstream infinity value 
of the tangential velocity. This result is also shown on 
Fig. 4; it appears that this is a better approximation to 
use than the downstream infinity values of the flow 
angles, especially at the root and tip. 

Velocity values at z = +0.17, are also shown, and it 
may be seen that these are about midway between disc 
and infinity values. This demonstrates that the decay 
rate of the perturbation components is rapid as the 
fluid moves away from the blades. 

Except in the expression for c,/Cm, the series for the 
perturbation velocity components converges rapidly. 
Six terms were actually calculated so that c,/c, and 
Co/Cm Might be accurately known. The values obtained 
by taking only two terms in the expansion are com- 
pared with the more accurate calculation on Fig. 4, and 
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Fic. 4. Comparison of exact with approximate solutions just 
downstream of the actuator disc. 


it may be seen that the added refinement is unnecessary 
+knz eT 


especially since the exponential e soon kills high. 
order terms away from the disc. The error involved in 
using Eq. (51) to calculate c,/c, cannot be more than 
halt the difference in the first series term—i.e., more 
than 0.002 in this case—since all the other series terms 
are virtually identical upstream and downstream. 


THE SECOND Exact SOLUTION 
Returning to the generalized equation for incom. 
pressible flow, Eq. (18), 
dH /dy = (1/r*)(0(d0/dy) + nr] 
suppose that nr is a function of 7 only. Differentiating 
with respect to s, 
d? FT Ow 1 d?[(1/2)0?] oy 
dy? oz r? dy? Oz 
Since y is a function of y and z, this equation can only 
be satisfied everywhere if 
d?HT/dy? = d*(0")/dy? = 0 


Integrating and again using nondimensional con 


stants, 
dH acy, H ™ 
= roe = on a a oo a3 i 
dy ry Cm” Cul i* 
d6 76a. \* y . 
7] = 12m ee ( = 2a. 2 + ay (33 
dy VrrOm Cm! 


Substituting the above into Eqs. (18) and (19) 
Woz + Vrr — (l/r), = gr = Cy(aiR? — az) (54 


Suppose, as before, that 
Y= bot Dery,’ (23 


Then 


Cm(Q1R? — a2) = Po, — (1/r)Wor + 


+knz 9 , 4 ' 

> e [LRnWn = Raa” =" i War 

n 
For this to be satisfied everywhere, 

9 , , l / 
Ra Ws a V nr <— WV nr = 0 
Z 

Putting A, = k,7r,, the solution for y’ has exactly the 


same form as in the case previously considered—Eq. 
(31). Also, 

Cmlai(R)? -- a>] = y.,, — (1/r)¥ 
may be solved to give 
ado 


=" Ri“ R?(210g,R— 1) —— R? (55 


T:"*Cm 8 4 Z 


< 


ds is another constant of integration, but it is not inde 
pendent of the others, being fixed by Eq. (26). Eqs. 
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FLOW THROUGH ANNULAR 


(93), (31), (52), (53), and (55) comprise full solutions for 


the flow. 
By differentiation, using the relations (12) and (13), 


Cc e/ . ~p 
~= > = bd,e7" Bi(AgR) (56) 
Cm n 
P ay - + n2/1 
&“__—-R?+a,+alog,R+ de "bh, Bo(\,R) 
9 
bs ~ n 
(57) 
The vorticity components are then given by Eq. (54) 
and 
r& = —A2(C,Cm/C9) (98) 
re = —A2(CCm/ Co) (59) 


We now have a second full solution, with, as before, 
three constants to determine the velocity distribution 
at infinity plus an additional constant aj, since varia- 
tion in stagnation conditions is permitted. The per- 
turbation effects may be treated in exactly the same 
way as in the first case, the values of \ depending only 
on the hub: tip ratio and the values of } being fixed by 
matching at the actuator disc. 

In matching two flows of this kind, the axial velocities 
at the dise will be exactly halfway between the in- 
finity values at the same radii for reasons similar to 
those given in discussing the case of uw, = wa¢ in the last 
section. This might be explained physically by con- 
sidering the flow to be caused by sheets of vortices 
shed by the blades that constitute the disc. Only the 
tangential component or ring component of these vor- 
tices will affect the axial velocity at the disc, and, since 
the strength of the vortex ring component at a given 
radius is the same everywhere, the net effect at the disc 
must be exactly half the effect at infinity. 

For constant stagnation enthalpy throughout the 
flow, a; = 0, and thus Eq. (54) reduces to 


nY = —A2xlm 


which is constant. This is a limiting case of the JJ = 
constant flows mentioned in the first section of the 


paper. 
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Types of flow with rc, = constant are obtained when 
dz = 0; then Eq. (54) becomes 


n/r = ACm/?r? 
which is constant, and 
C:./Cm = As — (1/2)a,R? (60) 


representing a limiting case of the flows with rc, = con- 
stant and »/r constant along a stream line discussed 
earlier. ‘ 

Although the tangential velocity distribution is the 
same as free vortex in this flow, the axial velocity is 
not the same because of the variation in stagnation en- 
thalpy with radius. The vorticity at a given radius is 
constant, and the flow consists of trailing ‘“‘smoke ring”’ 
vortices traveling along the stream surfaces, but the 
stream lines do not follow the vortex lines. 

From Eqs. (52) and (53), 


IT ( “itt) ay ( rCo ) cece 
= la; — 4 (61) 

og 2a2 2dae\P Cm 
Now the change in // along a stream line passing through 
a rotor is proportional to the change in rc. Since 
7¢o does not enter linearly into Eq. (61), it is not pos- 
sible to use this equation for flows both upstream and 
downstream of a rotor. Hence, for flow through a 
rotor, these solutions are restricted to cases for which 

H = constant. 

No numerical examples are shown for this second 


exact solution. 
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Remarks on the Interaction Between Parallel 
Streams of Subsonic and Supersonic 
Velocities 


Mac C. Adams 
Langley Aeronautical Laboratory, N.A.C.A., Langley Air Force Base, 
Va. 


January 3, 1950 


T° A RECENT PAPER, Tsien and Finston! presented an elegant 
treatment of the interaction between parallel supersonic and 
subsonic streams. 

For the case of the inclined wall in the subsonic region, there 
appears to be a possible error in the calculation of the pressure at 
the interface, Eqs. (54) and (55) of the reference. 

One might obtain the results for the inclined wall in the follow- 
ing manner: 

Using the notation of reference 1, the subsonic potential for the 
configuration shown in Fig. 1 can be written 


6 U2 1 ‘ y) 
"= ‘ V1 7; Pie sinh A Vl—- M,*(y of b) x 
: dd ' = 
(sin Ab cos d x) 2 + U; cosh AV1 — M2?(y + b) X 
6 


[a2(A) sin Ax + be(A) cos Ax]dA (1) 
The inclination of the flow at the wall is then given by 


2 6 WE hes d 
= sin Aob cos Ax 
obr 0 


“ 


~ 


6 
= [1(x + ob) — 1(x — ob)] 
2ab 


Now, if o — 0, Eq. 2 will give the unit impulse function—i.e., 
the limiting process gives a unit source on the boundary with a 
vertical velocity on the boundary of infinity at x = 0 and zero 





J 


M,>| -ob ob 





M,<| T ” 











T 


Fic. 1. Configuration considered in Eq. (1). 


elsewhere. The condition corresponds to a step on the Wall 
Certainly such a boundary condition is too severe for the limita. 
tions of the linearized theory; however, distribution of these steps 
would lead to solutions for an arbitrary boundary shape. For 
the case of the inclined wall, a uniform distribution is necessary 

Proceeding in a manner identical with that of reference 1, th, 
solution for the stepped boundary gives, for the pressure at th 
interface, 


(Ap): * 6 cos (0/2) oe §(O/x 
(1/JaUr ea — M2) (M2 — 1) coshé 


Now, integrating with respect to x leads to the solution for the 
inclined wall, where the quantity 6 is replaced by the wall angle, 
This results in 


( Ap),* 2e - 
: =~ - " 1 — 2cos — F3(é; 0) |], €>0 4 
(1/2)p. U1? VM? — 1 2 
(Ap)i* te ae 
= r cos ~ F,(é; 0), —& <0 
» 


(1 2) p, Ui? V M2 — | 2 


One might further suspect an error in reference 1 from the fact 
that negative pressures occur at certain points along the inter. 
face, Fig. 10 of the reference. 
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The Extension of the Prandtl-Glauert Rule 


E. V. Laitone 
Associate Professor, University of California, Berkeley, Calif. 
January 25, 1950 


—— APPLICATION OF THE Prandtl-Glauert Rule to three- 
dimensional flow and wind-tunnel wall effects has finally 
reached a satisfactory agreement in each case. An examination 
of most of the work prior to 1946!~3 indicates the inherent diff 
culties, mainly in not satisfying the proper boundary conditions 
and erroneous conclusions in applying the Prandtl-Glauert Rule 
By carefully considering the boundary conditions, the correct 
compressibility correction factors were finally obtained by 
Gothert and others.*~*® I believe the simplest rigorous deriva 
tion is the one I gave in references 7 and 8. _ In addition, this 
method permits a direct evaluation of the actual velocity distribu- 
tions on or about slender bodies. There has been some criti- 
cism of this proof because the perturbation potential was in the 
form described as Method I by Goldstein and Young.' How- 
ever, the correct boundary conditions were used in each case s0 
the same correct answer would be obtained independent of the 
choice of any constant factor multiplying the perturbation po- 
tential. For clarity and completeness, however, I should like to 
present the following as being the simplest rigorous approach: 


(1) Bopy oF REVOLUTION 


The linearized potential equation for compressible flow with 


axial symmetry is 
B*brz + drr + (1/r)b- = O | 
where for subsonic flow 
f 
B=V1— M 


and the fundamental potential for a source of strength F located 
at the point ¢ on the x axis and satisfying Eq. (1) is 
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o(x,r) = —(F/4n) }1 [(x — £)? + (6r)?] ) 
Then, fora continuous distribution of sources along the x axis from 
+= (tox = L, we have 
l I f(&) dé 
d(xn=- : : : (2) 
tr Jy [(x — £)? + (Br)?| 


The boundary condition to be satisfied at the body surface is 


fo = or (l + gz) = (or a O(drx) (3) 
Eqs. (2) and (3) can be combined as 
> a] 1. a 
Br , T(E) dé f(x) 
r ¢ — ; =~ ea 
tr oe WS = Ss)? > Ue) 2rr 
” r 
B2f"(x) rin + O(r) (4) 
L 
where the expansion of the integral expression is obtained by 
means of a Taylor series, as shown in reference 7. Therefore, we 
have 
fix) = 2arrory’U + O(ro4 In ry) = US'(x) + O(ro4 In 79) = (5) 
where S = mrp? is the cross-sectional area of the body of revolu- 


tion and the remainder term is exact only in the limit as 6 ap- 
proaches unity. It is seen that, at least to the first order, the 
required source strength for a given body shape is independent of 
compressibility effects, but the resulting velocities in compres- 
(2) and (5 


are greater, as shown by combining Eqs. 


the total velocity component in the direction of the 


sible flow 
For example, 
vaxis, for any point on or about the body, is 
1 L(x — &)S'(t) dé) 
[(x — )? + (Br)2}°/? f 


(6) 


u(x,7r) = U+ ¢ = ud4 ee 


which is the same correct result as given in reference 7. 
It is seen from Eq. (6) that the first-order wind-tunnel wall 

correction in compressible flow would vary as 1/8* for a body of 

revolution, since, for r >>(x — £), 

(x — £)S’(&) dé ( 


teri? J, 


As shown in reference 7, the surface pressure coefficient on the 


body of revolution is given by 


ste os ; 
Cx, ro(x)] = —2 "Tia (ro)? + O(ro4 In ro) = 
1 f% (& — x)S'(é) dé 
‘ ; 2 : : 7 — (ro)? + O(rot In ro) (8) 
2m 0 [(é — x)* + (Bro)? | 


the remainder term is exact only in the limiting 
case as 8 approaches unity. In this case, corresponding to in- 
compressible flow, Eq. (8) provides a satisfactory solution for the 
pressure distribution on either a sharp or founded nose body of 
revolution, provided the thickness ratio is less than 10 per cent. 
As shown in reference 7, we may expand Eq. (8) by means of a 
Taylor series to show that at any point on the surface of a sym- 
metrical body of revolution the first-order compressibility effects 


where, again, 


are given by 


Cry = Cp, — (S”"(x)/m) In B (9) 


for a body of 
as the selected 


This shows that the compressibility correction 
revolution is affected by the body shape as well 
location on a given body. 


(2) Two-DIMENSIONAL PROFILE 


A similar analysis is easily made for two-dimensional flow.® 
In this case the linearized potential equation for compressible 


flow is 


B°*br2 + dyy = O (10) 
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for a continuous source distribution along the x axis from 
x = Otox = C, the potential satisfying Eq. (10) 
1 c 
o(x, y) = f(é) In V (x — £)2 + (By)? dé (11) 
2a fi 
0 
and the boundary condition on the profile surface is 
yo’ = dy/(U + dz) = (by/U) + O(bydr) (12) 
Therefore 
, Bey : f(E) dé _— 
Uy’ = by = = * = Bf (x) + O(y?) 
2r Jo (x — &)? + (By)? 2 
so that 
f(x) = (2U/B)y0'(x) + O(y?) (13 


It is seen that now the source strength used for incompressible 
flow must be increased by 1/8 in order to represent the same pro- 
file in a compressible flow. 

Combining Eqs. (11) and (13), we derive the expression for the 
total velocity component in the x direction for any point on or 
about the body as 


ad Gi 


I (x — "(£) dé 


( py 3 


, E)y 
u(x, y) = UJ 1+ = (14) 
7B (x — £)? + 


Again, the first-order wind-tunnel wall compressibility correction 


would vary as 1/8*, since, for y >(x — &), 


( 


1+ (x — E)yo'(E) dé (15) 


u=U ses 
ry*p* 0 


However, on the body or along the stream line leading to the 
stagnation point, the compressibility correction varies as 1/8, 
since, for y ~ 0, 
od 4 4/> » 
© yo’ (E) dé 


1+ (16) 


= 
ll 
~ 


rB 0 z-€ 


It must be noted that for u values on the profile the denominator 
vanishes when & = x, so the Cauchy principal value of the in- 
tegral must be used. In this case, Eq. (16) provides a satisfac- 
tory surface pressure distribution for sharp nose profiles less than 


mek hs 


G=- = lim 
For rounded-nose profiles, it is necessary to use Eq. 


10 per cent thick as 
o (E)dE 

: (17) 
l rp . 


(14) for a 


satisfactory solution. 
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Stability of Alclad Plates 


P. P. Bijlaard 
Professor, Institute of Technology, Delft, Holland, on leave; 
Associate Professor, Cornell University, Ithaca, N.Y. 


January 21, 1950 


T° REFERENCE 1, the stability of Alclad plates, consisting of a 
core of high-strength aluminum alloy which is covered on each 
side with a coating of almost pure aluminum of very low strength, 
is investigated. The differential equation for the plates is derived, 
starting from stress-strain relations given by Ilyushin. The latter 
relations were derived from the present writer’s exact theory? by 
making the assumption that for the elastic deformation Poisson’s 
ratio is 0.5. Although this assumption does not cause much dis- 
crepancy with pure compression, it proves to give discrepancies 
from 20 to 30 per cent in the plastic reduction factor » in case 
of biaxial compression. This follows directly by inserting in the 
present writer’s stress-strain relations (21) or (2la) in reference 
2, with B = po/p, = 1, a value »y = 0.5 instead of the actual value 
of about 0.3. 

Moreover, it is easy to write down directly the differential 
equation of Alclad plates from the differential equation for homo- 
geneous plates,? 


O'w Ow O'w 
El\ A + 2(B + 2F) + D 


Ox! Ox*dy? oy! 
V O*w * V 21 0 (1 
l Ny - = ) 
* Ox? > dy? 


where the axes are assumed in the directions of the principal 
stresses p; and pz by which the plate is loaded, while Ny = hp, 
and Ny = hpe. 

The last two terms to the left in Eq. (1) refer to the deflecting 
forces exerted on a plate element by the compression forces Nx 
and Ny and consequently they remain the same for the Alclad 
plate. The first term to the left, however, refers to the restrain- 
ing force given by the transverse shearing forces. In the Alclad 
plate this term will consist of the sum of the pertinent terms for 
core and coating. With a core thickness / and a thickness ah 
of each coating, the moments of inertia of core and coating are 
h§/12 and [(1 + 2a)* — 1]h3/12 = gh3/12, respectively. Hence, 


the differential equation of the Alclad plate becomes 


B Ow n Otw O'w 
El, Ass 7 + 2(Ba + 2Fa) sae - 
x 


Ox*dy? : oy! 
N O*w -N 20 \@ 
d Ny = ( 2) 
* ox? * oy? 


in which 


Ig = h'/12, Aa = Ac + gdf 
Ba = Be + gBs, Da = De + gDy 
Fa = Fe + gfy 


The subscripts a, c, and f refer to Alclad plate, core, and coat- 
ing (faces), respectively. In the elastic domain 


A, = Da = Ba + 2F, = (1 + g)/(1 — v?) 


For a plate that is compressed in longitudinal direction and 
simply supported at the unloaded sides (Case 2 in reference 2), 
where for a homogeneous plate the plasticity reduction factor is 


(reference 2, Table 5): 
‘AD + B+ 2F ; 
ge = = 0.455 (\/AD + B + 2F) 


OF 2, (1 — pv?) 
it follows immediately that, for the Alclad plate, 


/ ~— 
VAD. B, 2F, 0.455 P ‘ z 
. T Be T SPs (V4qD, + Ba + 2Fa) (3) 
2(1 + g)/(1 — v2) l+g 


In the same way the formulas for 7 = og/og in Table 5 of 
reference 2 for Cases 1, 4, 13, 6, and 9 are valid for Alclad Plates 
if all values A, B, D, and F are given the subscript q and al 


values og/oz are divided by (1 + g). 
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Another Note on the Velocity of Sound 


C. B. Millikan 
Director, Guggenheim Aeronautical Laboratory, California Institute of 
Technology, Pasadena, Calif. 


January 14, 1950 


Is A RECENT NOTE,! Morduchow attempts to show that “the 

well-known differential expression a? = dp/dp for the velocity 
of sound cannot, as is sometimes attempted, be rigorously derived 
from the steady one-dimensional differential equations of non 
viscous flow.’”’ Unfortunately, his reasoning is faulty and his 
erroneous conclusion is likely to confuse unnecessarily the many 
students who have correctly derived this expression by the pro- 
cedure he believes to be incorrect. 

The three basic equations of the problem may be written, ina 
form slightly different from that used by Morduchow, as follows: 


du* = —(2 p)u*dp 
du? = —(2/p)dp 
du? = —([2y/(y — 1) ]d(p/p) 


where the gas is assumed to follow the perfect gas law. If we 
can find the conditions under which these equations are consistent, 
then these conditions will correspond to a one-dimensional steady 
flow in which uw, p, and p have constant values on cither side of 
a discontinuity across which there are infinitesimal discontinuities 
du, dp, and dp. The ratios dp/dp, du/dp, etc., are meaningless 
(mathematically 0/0) except at the discontinuity where they are 
well defined. 

For the second and third equations to be compatible we must 
have dp/p = [y/(y — 1) ]d(p/p), from which immediately fol 
lows p/p’ = const. This tells us that the gas flows isentropically 
across the discontinuity; thus, the entire flow is isentropic 
For the first and second equations to be consistent we must have 
u* = dp/dp, and combining this with the condition of isentropy, 
we find u? = y(p/p). 

Thus, the condition for the existence of a stationary infinitesi- 


mal discontinuity in a one-dimensional flow -with velocity « is 
u? = dp/dp = y(p/p) 


and this velocity we define as the velocity of sound, a. Mordu 
chow is led to what he calls a contradiction by taking the differen- 


tial of this relation, obtaining du? = yd(p/p). This is obviously 
inconsistent with the third of our original equations (unless y = 
—1). But the relation u2 = y(p/p) is not a general functional 


relation between u, p, and p which holds for all values of these 
variables. It merely gives a specific relation between them which 
must be satisfied for a stationary discontinuity to be possibl 
Accordingly, the relation cannot be differentiated and the expres 
sion du® = yd(p/p) is meaningless and incorrect. The simple 
“‘steady-motion”’ derivation of a? = dp/dp, which has seen such 


yeomen service, seems to be respectable after all 
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On Downwash Corrections in Unsteady Flow* 


John W. Miles 


Department of Engineering, University of California at Los Angeles 
December 16, 1949 


Seema HAS STATED that the most important unsteady flow 
effect in the calculation of dynamic stability derivatives 


generally may be identified as the ‘‘complex downwash lag’ 
factor, Which will enter in the determination of wing-tail inter- 
ference. While we have taken exception to the generality of his 
conclusions,? it does appear that the downwash lag will be im- 
portant, and we wish to advance here a simple method for its 
calculation 

Assuming the harmonic time deepndence exp (dwt), the ve 
locity potential (@) at the upper and lower surfaces of the vortex 
sheet behind a plane wing must exhibit a dependence on the 
streamwise coordinate (x, positive downstream) which will ensure 
continuity of the pressure. If we assume that the vortex sheet 
remains in the plane of the wing (z = 0), this requirement is satis- 


. 


fied by taking 


0 ) I ee a) | 
x, y,0 +) = +¢,(y) exp 4 — > 
Q(x y pi o CX] } 1 b { 


k = wh/U (2) 
where U is the freestream velocity, } is a reference length, y is 
the spanwise coordinate, x\(y) is the coordinate of the trailing 
edge, (x, y, z) is the velocity potential, and ¢)(y) is the value of 
this potential at the trailing edge of the upper surface. Calcu- 
lating the pressure at the wing and evaluating the section lift and 
moment (about x = 0) coefficients (c; and cm), we find 


ike(y)em(y) . 
(3) 
b 


cand ¢m are referred to c and c?, respectively, where c is the chord. 
Now, if c,,(y) represents the lift coefficient calculated in the ab- 
sence of unsteady flow effects and ¢o(y, 0 +) the corresponding 
velocity potential at the vortex sheet, we may write (assuming 


1 
di(y) exp [kx (y)/b] = Ul ey) - 


kto be small) 
cly) — (ike(y)em(y)/b) = Co(y) exp [ik6(y)] (4) 
o(x, ¥, 0 +) = go (y, 0 +) exp {ik[d(y) — (x/b)]} (5) 
Having the potential ¢(x, y, 0 +), the calculation of $(x, y, 2) is 
reduced to a classical boundary value problem. In principle, 
the solution of this problem is straightforward, but, in practice, 
it is likely te prove cumbersome; moreover, the introduction of 
such important factors as viscosity and warping of the vortex 
sheet leads to even more serious complications. However, if 
& (x, y, 2) represents the downwash field in the absence of un- 
steady flow effects but including experimentally or empirically 
determined corrections, it appears reasonable to assume that the 
downwash field in unsteady flow, at locations not too remote from 

the vortex sheet, may be approximated by 

: ada! ; ) . 
e(x, y, 3) = €o(x, y, 2) exp { ik [6(y) — (x/b)]} (6) 
Consider now the application of the result (6) to a pitching os- 
cillation of amplitude ap about the axis x = ab. The local angle 

of flow at the tail then may be written 


iy a 


aett. = aol l + ik[(x/b) — a}} — €7 (7) 





where er is the wing downwash at the tail. Substituting e from 


Eq. (6) and differentiating with respect to ao, we have 


dary, ' a) + ik (: 1 dn) (% ) a 
- = = 5 Sei 
dex day ; \ dao b : 
(a — 5) (“)] + O(k?) (8) 
dag 


* This work was carried out at the Aerophysics Laboratory, North Ameri 


can Aviation, Inc., during the summer of 1949 
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where e7, is implicitly dependent on the coordinates of a point 
on the tail and 6 on y. Now, in the absence of downwash inter- 


ference, the term ap in ae, would lead to a moment derivative 


Cutq,+ ikCuy, the latter term being due to unsteady flow effects, 
while the term skay|(x/b) — a] would lead to the derivative 
Cg. Hence, in the presence of interference, we write 
dC, ; dey ) a + ik , dey . 
= - Cm iR Ca. 4 
day da - daw : 


dery\ ? der, : 
1 — c M é + (a — 5) ¢ Vo +- ()(R?) 9g 
day dav . 


where e7, and 6 now are assumed to have been appropriately av- 
eraged over the tail. A similar result holds for the lift deriva- 
tive, and entirely analogous results may be obtained for the de 
We remark that the cor 
applied to Cu, 


rivatives with respect to angle of yaw. 


rections [1 — (der,/dag)] and [1 + (de7r,/dag)] 


and Cy, respectively, are standard, but the introduction of the 
term Cu;, together with its downwash correction and the term 
(a — 5)Cmuy, generally is neglected. The last term is presum- 
ably the ‘‘downwash lag”’ term referred to by Statler. It would 
appear that the usual analysis of downwash in unsteady flow* 
implicitly assumes 6 = a. 

In applying the result (9), the derivative (de7,/day) may be 
calculated from steady flow theory or measured in the wind 
tunnel, so that it only remains to calculate an appropriately aver- 
If we take this average as 

1/2 


aged value of 6(y). 


] 
S J _us 


c(y)6(y)dy (10) 


where 5S, c, and / are the area, chord, and span of the wing, and 
substitute 6(y) from Eq. (4), retaining only the first power of k, we 


have 
9 r 
1 ial fcigly) + cig ly) + [e(y)ema,(y)/b) | 
§=- c(y): - dy (11) 
R —|/2 ) Clay) f 
Assuming Cig, to be constant across the wing and taking x = 0 
at the aerodynamic center of the wing, we obtain 
§ = [(Crg + Cra)/Crq, |wing (12) 


For conventional wing-tail configurations, 6 will be negative. 
While it would have been a better approximation to have carried 
out the integration of Eq. (10) over the tail rather than the wing, 
it appears that the result (12) should be reasonably accurate in 
those cases where 6(y) does not exhibit large variations across 
the wingspan. Thus, the results should not be applied to heavily 
swept wings without further investigation 
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u v = (1) 
Ox oy pox 


may be interpreted as a balance between inertia and pressure 
forces. 
The Navier-Stokes equation 
Ou Ou 


+ + pV 2u (2) 
Ox oy p Ox 


is a balance between inertia, pressure, and viscous forces. 
Prandtl’s boundary layer theory assumes that 


viscous effect 
inertia effect 


or that the Reynolds Number is extremely large. This leads to 
the familiar boundary-layer solution and holds for extremely small 
values of viscosity provided the length parameter in the Reynolds 
Number is also not small. 

Equations such as Eq. (2) have been the subject of some recent 
study.' It has been shown that, if the parameter that multiplies 
the highest term of the differential equation approaches zero 
(i.e., if y + 0 above), then the solution to the equation approaches 
the solution corresponding to the equation for » = O [Eq. (1) 
above |, except for a narrow region within which there is a rapid 
transition of—in this problem—velocity. 

The phenomena are dependent upon the number of boundary 
conditions required for a unique solution of a partial differential 
Thus, Eq. (2) being of a higher order than Eq. (1), it 


equation. 
the no-slip condition. 


requires an additional boundary condition 
For v — 0, one is tempted to neglect the last term of Eq. (2), but, 
if this were done, the order of the equation would be lowered, a 
boundary condition would be lost, and the solution (at least in 
the neighborhood of the boundary) would differ from the no-slip 
solution. 

“True” boundary-layer solutions 
the form of the differential equation is similar to Eq. (2) 
been discussed. Taylor and Maccoll,” in their analysis of the 
shock front of supersonic flows, obtain a boundary-layer equation 
that leads to the characteristic thin layer (shock thickness) as 
viscosity and heat conduction approach zero. The laminar flow 
heat convection problem? also leads to an equation that is identi- 
cal in form to the viscous flow equation. 

_ Friedrichs‘ discusses a differential equation that occurs in the 
theory of plates and shells and which physically may be thought 


of as a balance of 


that is, solutions in which 
have 


plate effect plus membrane effect = transverse load effect 
and considers the case in which 


plate effect 0 
membrane effect 


It is found that for this case the boundary layer is in the slope 
of the structure at the supports and, in effect, that the structure 
acts like a membrane everywhere except at the supports, at 
which place a boundary condition is lost and a rapid transition 
in slope takes place in a narrow region, this being the boundary 
layer. 

Another similar, although apparently unrelated, boundary- 
layer type effect occurs in the theory of plates.5 If we consider 
a long plate, simply supported along its long sides, free along the 
short sides, and uniformly loaded transversely over its area, it 
is found that at the short edges the plate ‘‘curls’’ over in a narrow 
region—the boundary layer for this case. The phenomenon is 
undoubtedly due to a boundary condition effect; some distance 
away from the short edges a condition of plane strain exists, 
whereas at the edges this condition must be given up. The 
boundary-layer condition described by Friedrichs may be 
thought of as one in which a condition of plane stress prevails 
throughout the plate but is given up at the edges. 
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Bromberg and Stoker® discuss a boundary-layer effect th 
occurs in the theory of shells which is due to the inclusion jy 1, 
plate theory of certain nonlinear, higher-order terms of the str;j, 
tensor. 

Two additional types of boundary-layer effects which are dy 
to singularities in the boundary conditions and which lead 4, 
characteristic thin layer regions will be described 


The first occurs in the Lame-Clapeyron theory of axially syp. 





metric deformations of elastic bodies. Consider a rotating circy. 


lar wheel. It is found’ that the tangential and radial stress djs 
tributions for the case in which there is a small central hole jn thy 
wheel are the same as the solutions corresponding to no hole 
everywhere except for a narrow region adjacent to the hole ; 
which there is a rapid change in the stresses. This is a boundary. 
layer type effect and is connected with the following singularity 
in the boundary condition; for the case of no hole, the boundary 
condition is: (a) The radial deformation is zero at the center 
of the wheel. 
(b) The radial stress at the edge of the hole is zero, Thy 


For the case of a small hole, the boundary condj- 
tion is: 
differential equation in deflection is the same for both cases, but 
in going from small hole to no hole there is a singularity in the 
boundary condition. The solution is similar in many respects 
to that of the small circular hole in an infinite plate subjected to 
tension stress. Also, there appears to be some connection with 
the St. Venant Principle. 

Another boundary-layer effect occurs in the direct compression- 


buckling loading of a rod. If we consider the ratio of effects 


buckling (bending) stress 
>0 


direct compression stress 


it is found that the governing parameters are the slope (6) and 
the slenderness ratio (//p) of the bar. The slenderness ratio being 
fixed in any given problem, it becomes necessary that the slope 
approach zero. 


The problem can be set up as follows: for zero slope, 


o = P/A 3 
for small slope, 
§o = (P/A) * [(Py)c/T] fa 
1 El(de/dS) = Py 4b 
in which 
= stress 
P = load 
A = area 
y = deflection 
c = distance from neutral axis of bar 


EI = stiffness 

dS = differential length along bar 

For extremely small 6, the stress (4a) is nearly the same as the 
stress (3) for all values of P less than Pguter. As P approaches 
and then exceeds Pguter, it is found’ that there is a rapid increase 
in o corresponding to a small increase in P. Thus, the boundary- 
layer ‘‘region”’ in this case is the load P, and the term that corte- 
sponds to the velocity of the viscous boundary layer is the stress. 
For small 6, a boundary condition on the slope at the top of the 
bar is required, whereas for zero slope no boundary condition is 
needed. Thus, in going from small slope to zero slope a boundary 
condition is lost. 

A familiar example of some interest—in that it indicates the 
importance of losing a boundary condition (dropping the highest 
order term of the differential equation)—in the boundary-layer 
type solution is the comparison of Stokes and Oseen flows.’ The 
Stokes flow is a balance of viscous and pressure forces and leads to 


vV *u = (1/p)(0p/dx) (5 


The Oseen improvement introdtces a perturbation inertia term 
(first order in the differential equation) 
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April, 1947 fuselage length is 34.5 ft. 
(B) Wings 


The span is 40 ft.; area, 240 sq.ft.; mean aerody- 
namic chord (M.A.C.), 6.75 ft.; and aspect ratio, 6. 


Appendix 


SPECIFICATIONS OF AIRPLANE USED IN EXAMPLE 
(A) Airplane—General (C) Horizontal Stabilizer 


The airplane is a low-wing, fighter type, powered by The span is 16 ft.; area, 45 sq.ft.; aspect ratio, 6; 
a jet-propulsion engine exhausting through the extreme elevator hinge centerline, at 75 per cent chord; and 
aft end of the fuselage. The gross weight is 10,000 Ibs., tail length (distance from c/4 of wing M.A.C. to ¢/4 
and, with a center of gravity location of 27 per cent of tail M.A.C.), 15 ft. 
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